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INTRODUCTION 


It  has  been  observed  empirically  (Gringorten,  2966,  1968)  that  fluc¬ 
tuations  in  upper-level  wind  velocity,  rainfall,  temperature,  and  other 
meteorological  variates  may  be  modeled  by  a  type  of  random  Markov  process. 
A  good  fit  is  obtained  over  a  broad  range  of  values  when  the  so-called 
Omstein-Uhlenbeck  process  described  below  is  employed.  Acceptance  of 
the  validity  of  such  a  model  permits  utilization  of  its  theoretical  pre¬ 
dictions  for  long-range  operational  planning,  with  on-line  predictions 
available  from  software.  Such  software  is  in  development  at  the  U.S.A.F. 
Geophysical  Research  Laboratory. 

The  0-U  process  is  the  simplest  random  process  having  a  stationary 
distribution  whose  sample  paths  are  continuous  functions  of  time.  It  is 
correspondingly  easiest  to  work  with.  Certain  passage  time  distributions 
needed  for  meteorological  planning  have  been  difficult,  however,  to 
calculate,  and  have  only  recently  become  available  (Keilson  and  Ross, 
1975). 


The  stationary  Gaussian  Markov  process  or  0-U  process  X(t)  has  the 
following  characteristics 

(a)  E[X(t) ]  =  0 

(b)  Var  X(t)  =  1 

(c)  cov[X(t),  X (t  +  t)]  =  e"l  I 

(d)  The  probability  density  function  is  standard  normal,  i.e., 


fX(t)(x) 


(e)  The  transition  probability  cumulative  distribution  function 


is 
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x  exp{ - (y  -  xfte_t)2/2(l  -  e'2t)} 

P[X(t)  s  x | X(0)  =  x  ]  =  /  - - - -  dy 

/2ir(l  -  e”2t) 

Many  related  variates  are  of  meteorological  interest.  Of  central 
interest  to  this  study  has  been  the  maximum  value  of  X(t)  attained  over 
an  interval  of  specified  duration  6,  i.e., 

M(e)  =  max  X(t) 
t<t'<t+0 

The  distribution  of  M(9)  relates  directly  to  the  passage  time  dis¬ 
tributions  previously  obtained  (Keilson-Rcss,  1975).  The  relationship, 
though  direct  and  straightforward  in  principle,  has  required  extensive 
analytical  and  numerical  effort  to  convert  to  useful  form.  The  results 
needed  have  been  attained.  A  concise  presentation  of  the  results  was 
given  at  Banff,  Canada  in  October,  1979  at  the  Sixth  Conference  on 
Probability  and  Statistics  in  Atmospheric  Sciences.  The  paper  given 
there  is  Part  II  of  this  final  report. 

The  Banff  paper  of  Part  II  is  a  precis  of  a  longer  paper  documenting 
our  results  in  full  detail  with  complete  tables  and  graphs  given.  This 
longer  paper  is  Part  I  of  this  final  report. 

The  programs  needed  to  make  the  numerical  results  available  on-line 
are  available,  and  will  be  developed  for  USAF  use  at  the  Air  Force  Geo¬ 
physics  Laboratory  (AFSC)  at  Hanscom  Air  Force  Base  in  Massachusetts. 
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1 ,  Introduction  and  Summary 

The  Gaussian  Markov  process,  also  known  as  the  Ornstein-Uhlenbeck 
process,  has  been  widely  employed  to  describe  fluctuations  in  physical 
phenomena.  It  has  also  been  a  basic  theoretical  and  computational  tool 
for  statistics.  References  to  such  0-U  applications  have  been  presented 
in  [10],  where  0-U  passage  time  distributions  were  described  and  tabulated. 

An  important  potential  application  of  the  0-U  process  is  to  meteoro¬ 
logical  prediction  and  operational  planning.  I.  Gringorten  has  reported 
extensively  [8,9]  on  the  empirical  validity  of  the  0-U  process  as  a  model 
for  fluctuations  in  temperature  and  other  climatic  variates.  For  meteoro¬ 
logical  applications,  the  range  of  such  vaTiates  and  the  persistence  of 
high  or  low  values  of  such  variates  are  of  considerable  importance.  This 
report  is  concerned  with  the  maximum  of  such  vaTiates  over  a  specified 
time  interval  predicted  by  an  Ornstein-Uhlenbeck  model  for  the  fluctua¬ 
tions.  More  precisely,  we  are  interested  in  the  maximum  value  M(0)  of 
the  process  over  an  interval  of  length  0,  when  the  process  has  its  sta¬ 
tionary  distribution,  the  standard  normal  distribution,  at  the  beginning 
of  the  interval.  We  study,  therefore,  the  sequence  of  random  variables 

(1)  M(0)  «  max  X(t') 

Ost  *se 

where  X(t)  is  the  stationary  0-U  process  with  distribution 

x  - y2/2 

(2)  Fe(x)  «  P[X(t)  s  x]  -  /  * -  dy  ,  0  s  t  <  •  , 

—  /27 
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and  covariance  function 

(3)  rx(t)  =  E[X(t)X(t+t)]  =  e”T  . 

Physical  processes  of  interest  described  by  an  0-U  model  have  three 
parameters  of  interest:  .  the  mean  value  m,  the  rms  amplitude  A,  and  the 
time  scale.  A  process  V(u)  with  mean  m,  amplitude  A  and  correlation 
function  Py(t)  *  cov[V(u),  V(u*t)]/A^  «=  e  T/^  would  be  described  by 

(4)  V (u)  *  m  ♦  AX (u/T) 

where  X(t)  is  our  0-U  process.  We  note  that  both  the  process  X,  its 
states  x  and  time  parameter  t  are  dimensionless.  The  O-U  process  has 
zero  mean  and  standard  deviation  (amplitude)  equal  to  one. 

The  distribution  of  the  maximum  M(0)  for  interval  length  6  is  related 
to  the  passage  time  distributions  previously  described  [10].  The  calcu¬ 
lation  technique  employed  there  evaluated  zeros  and  residues  in  the  com¬ 
plex  plane  of  the  parabolic  cylinder  functions.  Such  zero  and  residue 
evaluation  techniques  again  feature  in  tabulating  the  distribution  of 
M(0) . 

The  extreme  values  of  M(0)  for  large  0  have  been  described  in 
indirect  form  by  Darling  and  Erdos  in  the  form  of  a  limit  theorem  [4]. 

The  appropriate  form  of  the  limit  theorem  is  presented  in  Section  3. 

Such  limit  theorem  information  has  been  of  limited  value  for  prediction 
and  planning  because  such  limit  theoretic  behavior  has  been  known  to  set 
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in  quite  slowly,  and  when  the  limit  theorem  was  useful  has  not  been  known. 
Correspondingly,  the  error  in  the  limiting  distribution  has  not  been  known 
for  any  value  of  6. 

The  distribution  of  M(0)  and  the  underlying  theory  supporting  its 
tabulation  are  given.  The  distribution  F0(y)  of  M(0)  and  the  associated 
density  f0(y)  are  tabulated  and  plotted  graphically.  Comparison  with 
the  limit  theorem  of  Darling  and  Erdos  is  discussed. 

A  basic  feature  of  the  behavior  of  M(0)  is  that  for  large  values 
of  6,  M(0)  loses  its  variability.  Formally,  one  has 

(4)  M(0)  ^  /21og6  , 

i.e.,  M(0)  converges  in  probability  (or  distribution)  to  the  constant 
/21og0 .  Such  behavior  is  familiar  in  extreme  value  theory,  accounts  of 
which  may  be  found  in  [7]  and  [1], 
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2.  Basic  Theory 

A.  The  distribution  of  the  maximum  and  its  density 

Let  X(t)  be  the  stationary  0-U  process  [3]  for  which  the  density 
of  X(t)  is  standard  normal,  i.e., 

»>  w*  --%r',2/2 

and  the  covariance  function  is 

(2)  rx(i)  *  cov [X(t) ,  X(t  4  t)]  *  e"T  . 

Let  M(8)  be  the  maximum  of  X(t)  over  an  interval  of  length  0,  and  let 
F0(y)  be  its  c.d.f.,  i.e., 

(3)  Fg(y)  *  <  y]  *  P[X(f)  <  y,  0  S  t'  S  6]  ,  -“>  <  y  <  “  . 

We  wish  to  evaluate  the  distribution  Fg(y)  and  its  p.d.f. 

(4)  fe(y)  -  F0(y)  ,  -•  <  y  <  •  . 

We  also  seek  to  quantify  the  asymptotic  behavior  of  Fg(y)  and  fg(y)  as  6  +  •, 
and  compare  our  results  with  the  limit  form  of  this  distribution  of  extreme 
values,  given  effectively  by  Darling  and  Erdos  [4], 

The  methods  employed  in  [10]  to  calculate  the  passage  time  distributions 
for  the  0-U  process  are  applicable  with  some  extension  to  the  calculation 
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of  the  distribution  of  M(6).  The  desired  variate  M(6)  is  related  to  the 
stationary  O-U  process,  and  this  feature  of  stationarity  must  be  dealt 
with. 

To  understand  the  relation  of  the  passage  time  problem  to  the  maximum 
problem,  a  discussion  in  the  context  of  a  simpler  birth-death  process  will 
be  helpful. 

Consider  a  birth-death  process  N(t)  assumed  to  be  stationary  on 
states  {0,1,2, ... ,K}.  Let  M(0)  *  max  N(t').  Then  if  T  is  the  time 

Ost'se 

to  go  from  state  r  to  state  n  and  if  er  is  the  probability  that  the  sta¬ 
tionary  N (t)  is  in  state  r,  one  has 

n-1  K 

(5)  P[M(6)  ^  n]  =  l  P[T  <  0]e  ♦  l  e 

r=0  r=n  r 

n-1 

■  1  -  7  e  P[T  >  6]  , 

r=0 


so  that 

n-1 

(6)  P[M(6)  <  n]  *  l  e  P[T  >  e]  . 

r*0  n 

-sT  •» 

Let  E[e  rnJ  *  o  (s)  *  /  e“STs  (T)dT  where  s  (x)  is  the  p.d.f.  of  T  . 

***  q  m  xn  rn 

Let  5  (t)  *  PfT  >  tJ  be  the  survival  function  of  T  .  Then  one  has 

„  rn  rn  rn 

/  e  ^^rn  >  T^T  *  “  orn(s)]/s.  Hence  from  (6)  we  have 


(7) 


*  cfl  n-1  1 1  -  o  (s)  I 

/  e'S0P[M(6)  <  n]de  *  7  ej - ^ - } 

0  r*  0  r  s 
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The  0-U  process  may  be  regarded  as  the  limit  of  a  sequence  of  birth-death 
processes  and  the  structure  of  equation  (7)  then  carries  over  directly  to 
give 


(8) 


/  e”s6P[M(6)  <  y]de  =  jf 
0  /2tT 


X  e'x2/2 


l-o  (s) 
xy  J 


dx 


where  M(6)  is  the  stationary  maximum  for  the  0-U  process  and  o  (s)  is  the 

xy 

Laplace  transform  of  the  passage  time  density  from  x  to  y.  It  has  been 
shown,  however,  in  [5]  that 


D  (~X)  -  y 

(9)  oxy(s)  =  expKx*  -  y')/4] 


where  Dv(z)  is  the  Weber  function  of  order  v  ([12],  p.  323,  et  seq.).  Equa¬ 
tions  (8)  and  (9)  are  the  basic  equations  for  our  calculations.  As  shown  in 
Appendix  1,  the  function  D_s(-y)  is  an  entire  function  of  s  for  all  y,  with 
zeros  only  on  the  negative  real  s  axis.  These  zeros  are  simple  and  give 
rise  to  simple  poles  for  the  expression  in  (9),  and  the  locations  of  the 
poles  and  residues  contributing  to  P[M(0)  <  y]  for  the  inversion  of  (8) 
are  found  as  in  [10].  It  is  shown  in  Appendix  1  that 


(10) 


J  * 


-x2/2 

vs,dx 


,-v2n  p.,.j(-y)  a  y  -  »„w 

GZ  "D-s(-yi  /2tT  nil  5  4  Vy) 


where  (y)  and  @n(y)  are  positive  and  0  <  Aj(y)  <  A^ (y) -  - .  •  From  (8) 
one  then  has  (Appendix  1) 
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-y  /2  »  8_(y)  -x  (y)e 


(11)  F0(y)  =  P[M(6)  <  y)  *  S— -  J  JL—  e 

/2t  n=l  nl/' 


This  may  be  rewritten  as 


r  -xnCy)e 

(12)  F  (y)  =  F  (y)  l  p  (y)e 

n=l  n 


y  2 

where  F0(y)  *  /  — e"x  ^2dx,  and 
J  -®  /JtT 


(13) 


Pn(y) 


e*y  l2//2^  gn(y) 

Frt(y) 


*n(y) 


B.  The  ergodic  exit  time  density  sE(x) 

The  integral  appearing  in  (10)  relates  to  a  simple  entity  of  probabilistic 

F 

interest,  sy(x),  the  "ergodic  exit  time  density"  from  (-«■  ,y)  [11].  This  is 
the  p.d.f.  of  the  time  to  depart  from  (-®,y)  if  one  starts  with  the  stationary 
distribution  truncated  to  (~®,y)  and  renormalized  to  unity.  Thus 
sy(T)  *  *  **  [X ('t  * )  5  y.  0i  t'  s  t  | X CO)  s  y] .  Its  Laplace  transform  is 

then  given  by 

2/o 

y  -x  /  2 

(14)  0*(s)  *  (F  (y)}-1  /  £— -0  (5)dx  . 

y  u  -«  xy 


It  is  known  that  Markov  diffusion  processes  are  reversible  in  time  and 
E 

that  sy  (t)  for  such  processes  is  completely  monotone  in  x  [11].  It  then 
follows  that  the  integral  on  the  right  hand  side  of  (8)  is  equal  to  F0(y)§E(x) 

p  ®  g  y 

where  §y(t)  *  /  sy(t')dt'  is  the  survival  function  of  the  ergodic  exit 
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time  and  is  also  completely  monotone  in  t.  Hence  pn(y)  must  be  positive, 
and  this  positivity  is  borne  out  by  the  computation.  We  note  from  (12) 

oo 

that  l  p  (y)  *  1.  This,  too,  emerges  from  the  computation  and  provides 
n=l  n 

an  accuracy  check. 

The  density  of  M(0),  f.(y)  may  then  be  obtained  from  (12)  via 
numerical  differentiation.  Alternately,  (Appendix  1)  one  finds  analy¬ 
tically  that 


(15) 


Mfe(y)] 


Vy> 

s  -  Xn(y) 


2 


where  L  [f_(y)]  *  /  exp(-s8)ffl(y)d6.  Hence,  from  (10)  and  (14),  we  have 
s  e  Q  e 

jy2  7  E  F 

(16)  f0(y)  *  v'H  F‘(y)  {s*(e)|s*(8)}  . 


We  note  from  (15)  that  lim  si  (fa(y)) 

s-~  5  6 


fp(y)  *  e”y  implies  that 


(17) 


I 


en(y) 


1  S  +  Vy'  /s 


Since  from  (10),  (14),  and  (17), 


tsISy(9)l 


Oy(S) 


e'yV2//2 7 


Vy>  1 

s  *  An(y)_ 


e‘y2/2/^  1 
Vy>  /s 


s  ♦  • 


The  state  x  -  0  is  a  special  state  for  the  O-U  process,  and  many  of 
the  descriptive  densities  of  the  process  take  a  simpler  form  when  this 
state  is  involved.  The  c.d.f.  F&(y)  and  the  density  f0(y)  also  take 
simpler  form  for  y  =  0,  and  the  closed  form  answers  obtained  are  struc¬ 
turally  informative.  Details  are  provided  in  Appendix  2. 

The  ergodic  exit  time  density  Sq(t)  from  the  interval  (-*,0)  is  given 
(Appendix  2)  by 

(19)  s£(0  -  |  e'T(l  -  e’2Y1/2  . 

This  may  be  rewritten  as 

(20)  s  (t)  -  l  p  ne"nT  , 

c  n  odd 


where  pn  *  pn(0)  (cf-  (3. 8))  is  the  value  of  our  basic  positive  weights 
pn(y)  at  y  *  0.  The  pn  are  shown  in  the  appendix  to  be  given  explicitly  by 

(21a)  p2.  -  p2.(0)  -  0  ; 
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One  then  has,  asymptotically,  for  the  slowly  convergent  sequence  p^^. 


(22)  p2k+1  ~  (wn) ' 


.  n 


We  note  from  (20)  that  the  ergodic  exit  time  Tq  has  all  moments.  We  also 
see  that  s^ft)  is  completely  monotone  in  t.  A  direct  integration  of  (19) 
gives 


5q(t)  *  |  arcsin(e-T)  , 


an  arcsin  law  for  T^.  It  then  follows  directly  from  (23)  nd  (3.7)  that 


(24)  F0 (0)  «  P[M(0)  <  0)  «  i  arcsin  e-6 


Hence,  one  has  the  asymptotic  behavior 


(25)  P[M(6)  <  0]  ~  -J-  e'6  ,  6  » 


i 
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3.  The  Asymptotic  Distribution  of  M(9) 

A.  The  Basic  Limit  Theorem 

It  may  be  inferred  (see  Appendix  3)  from  the  ideas  and  results  in 
Darling  and  Erdos  [4]  that  the  interval  maximum  M(6)  satisfies  the 
limit  theorem 

(1)  lim  P[c(0)(M(e)  -  ;(©)}  <  x]  «  exp(-e'x}  =  G(x) 

e-*» 

where 

(2a)  c(e)  *  /21og6  , 

and 

(2b)  ic«)  ■  c(.) .  isfjip  -  . 

B.  Validity  of  the  Limit  Distribution 

The  appearance  of  the  extreme  value  distribution  G(x)  in  (1)  is  to  be 
expected,  of  course,  from  the  more  basic  results  of  extreme  value  theory 
as  developed,  for  example,  by  Gnedenko  [6]  and  explored  statistically  by 
Gumbel  [7].  Indeed,  the  nature  of  the  stationary  0-U  process  is  such 
that  the  sequence  of  values  (X(ak),  k  =  0,1,2,...)  are  normally  distri¬ 
buted  and  asymptotically  independent  as  a  +  ».  This,  of  course,  is  true 
for  any  stationary  Gaussian  process  Y(t)  whose  covariance  function  rY(t) 
vanishes  for  large  t.  It  might  then  be  anticipated  that  the  extreme 
value  theory  developed  for  discrete  sequences  of  independent  random  vari¬ 
ables  is  still  valid  for  the  continuous  time  0-U  process.  The  Darling-Erdos 
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theory  provides  the  proof  that  this  is  so.  It  does  not,  however,  say  for 
what  values  of  0  and  y  the  approximation 

(4)  Ffl(y)  *  P{M(6)  <  yj  =  G(c(6)y  -  c(6);(6)) 

implicit  in  the  limit  theorem  of  Eq.  (1)  is  valid.  It  also  provides  no 
information  on  the  error  contained  in  the  approximation  (4) .  The  tabulated 
values  of  F„(x)  and  associated  graphs  are  therefore  essential  for  predic- 
tion  and  planning. 


C.  "Deterministic”  Behavior  of  M(6)  for  Large  6 

Since  c(0)  -*•«*>,  as  6  +  •,  and  c(6)  -  c(6)  4  0  as  0  -*■  *»,  it  follows  from 
(1)  that 


lim  P[ |M(0)  -  c(0) j  <  tj  *  0  ,  all  e  >  0  , 

e-*«o 


i.e.,  that 

(S)  M(9)  -  c(0)  5  0  ,  e  •  . 

This  implies  that  M(0)  converges  in  distribution  to  the  constant  c(6)  as 
0  -►  *,  i.e.,  that  M{0)  loses  its  uncertainty  asymptotically.  One  also 
expects  that  M(0)  converges  in  time  to  c(0)  for  every  sample  path  with 
probability  one,  i.e.,  that  the  strong  law  holds.  One  knows  that  M(0) 
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has  all  moments,  and  that  the  asymptotic  law  G(x)  for  c(6){M(6)  -  5(0)}  has 
all  moments.  From  this  and  the  numerical  information  on  fQ(y)  one  expects 
that  M(6)  will  have  all  moments  for  all  6.  A  proof  seems  difficult. 

Since  5(0)/c(0)  -*■  1,  as  6  -*■  •  ,  c(0)  describes  both  the  displacement  of 
M(0)  from  zero  and  the  scale  of  the  deviation  of  M(8)  from  its  mode  as 
measured,  e.g.,  by  the  half-width  of  M(0)  about  its  mode.  The  linear 
relation  between  the  mode  and  such  a  scale  factor  is  evident  from  the 
graph  of  the  density  ffi(y). 

A  simple  table  of  values  for  c(0)  and  5(0)  is  shown  in  Figure  1, 
demonstrating  the  extreme  slowness  with  which  c(0)  increases  with  0. 


C  (6) 

C(6) 

0 

1 

0.08 

1.65 

2 

1.89 

7.39 

3 

3.06 

9. OOxlO1 

4 

4.12 

2.98xl03 

5 

5.14 

2.68xl05 

6 

6.15 

6.57xl07 

FIGURE  1 


The  mean  and  variance  for  the  distribution  function  G(x)  may  be 
obtained  from  its  moment  generating  function  (cf.  [7]) 

•  SO 

Y(s)  *  /  e  sxdG(x)  *  /  esxe*xexp{-e"x)dx 

••  —SO 

so 

-  /  t”se_tdt  «  T(l-s)  . 

0 


(Sa) 
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We  then  have  *  -r'(l)  «  0.577215 _  (Euler's  constant),  and  from 

aG  =  (log  y(s)]q  we  find  o^,  =  n  / 6,  i.e.,  *  1.28255.  Equation  (1) 

may  be  employed  to  write  for  M(0)  the  approximate  distribution  valid  for 
large  0 


(5b) 


M (0)  »  C(0)  + 


2 

where  £G  has  mean  and  variance  wG<  and  oG  given  above.  We  find,  there¬ 
fore,  that  for  large  0 


(5c)  E(M(0)]  *  c(0)  ♦  wG/c{0)  ; 

°G 

(5d)  °M(0)  *  EM  * 

The  density  g(x)  of  SG  is  e  exp{-e  }  and  has  its  mode  at  0.  It 
then  follows  from  (5b)  that  for  large  0  the  density  f0(y)  of  M(0)  has  its 
mode  at  c(0) ,  i.e. , 


yM0DE(e)  *  c(e)  *  e  »  1  . 

It  is  of  interest  to  compare  the  asymptotic  behavior  of  M(0)  with  that 

of  M.  *  max  X.,  where  X.  are  i.i.d.  and  standard  normal.  In  that  case, 
lsjsk 

one  finds  that  converges  in  probability  to  *^21og  k  with  correction  terms 
of  a  character  similar  in  form  to  that  for  M(0).  Details  may  be  found  in 
Cramer  [2]  and  in  [1].  That  and  M(k)  both  go  to  /21og  k  in  probability 
is  in  keeping  with  the  fact  that  M(0)  for  the  0-U  process  has  a  "relaxation 
time"  equal  to  unity. 


1 
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D.  Relation  of  the  Asymptotic  Distribution  of  M(8)  to  X  (y)  and  Pn(y) 

For  every  value  of  y,  one  has  the  relation 

(6)  F0(y)  =  P[M(0)  <  yj  =  P[M(0)  <  y] P [T^  >  6] 
i.e. , 

(7)  F0(y)  =  FQ(y)Sy (6) 

E 

where  is  the  ergodic  exit  time  discussed  in  Section  2,  part  B.  This 

may  be  inferred  immediately  from  (2.8)  and  (2.14).  For  large  y, 

P [M (0)  <  y]  -*■  1  and  the  asymptotic  behavior  of  M(6)  is  contained  in 

-E 

the  survival  function  Sy(0)  of  the  ergodic  exit  time.  From  Equations 
(2.12)  and  (7)  we  have 

-E  %  *Vy)0 

(8)  S‘(fl)  =  l  p  (y)e 

y  n=l  n 

A  graph  of  the  eigenvalues  *n(y)  plotted  against  y  is  given  in  Figure  2. 
It  is  seen  that  all  the  eigenvalues  X^(y)  decrease  with  y  for  all  real  y, 
and  that  for  the  principal  eigenvalue  Xj(y),  one  has 

(9a)  lim  X.(y)  =  0  ;  X. (0)  *  1  . 

y-nx>  J 

For  the  other  eigenvalues,  one  finds 

lim  XR(y)  «  n-1  ;  An(0)  «  2n-l  , 


(9b) 
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and  one  notes  that  for  0  <  y  <  •»,  the  values  A^(y)  have  approximately 
linear  spacing  (Fig.  2). 

Ke  have  seen  in  Section  2  that  Pn(y)  >  0,  and  J  Pn(y)  =  1  for  all  y. 

It  is  also  found  numerically  that  p^(y),  the  mass  associated  with  the  maxi¬ 
mal  eigenvalue,  increases  monotonical ly  with  y  for  all  real  y  and  that 

(10a)  lim  Pj(y)  =  1  . 

y-Ko 

Correspondingly,  one  finds 

(10b)  lim  Pn(y)  =0  ,  n  i  2  . 

y-*oo 


A  graph  of  Pn(y)  obtained  from  the  computer  is  shown  in  Figure  3. 
run  or  against  r  ro*  ■>  -  K»«  rim  „  p  (y)  „  F0R  „  .  I(1)t 


FIGURE  2 


FIGURE  3 
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From  (7),  (8),  (9a)  and  (10a)  one  obtains  the  asymptotic  relation 

-x,(y)e 

(11)  Fe(y)  ~  Pi(y)F0(y)e  .  y  -*•  *  , 

and  the  corresponding  approximation 

-*.(y)e 

(12)  Fe(y)  s  Pl(y)F0(y)e  1 
It  is  found  that 

*Ai(y)e  4 

(13)  0  s  Ffi(y)  -  P1(y)Ffl(y)e  s  10  *  ,  6  >  12  ;  —  <  y  <  » 

i.e.,  we  have  accuracy  for  (12)  to  four  significant  figures  for  all  real 
y  when  6  >  12  or,  say,  c(6)  2  2.5. 

It  may  be  noted  that  since  (6)  is  true  for  all  6,  y,  we  may  replace 
6  by  p0yT  where  u0y  «  E[Ty].  Hence 

(14)  P(M(u0yt)  <  y]  «  F0(y)P[T^/U()y  >  t]  . 

But  E[Ty]  ~  E[T0y),  y  -*•  •  and  P[(Ty/nEy)  *  -*■  e‘T,  as  x  •  [4,  11]. 

Hence 

(15)  lim  P[M(t)  <  y^(y/t)]  =  e'T 

y**  y 

INV 

where  yQy  (x)  is  the  inverse  function  to  x  *  v0y,  which  is  increasing  on 


(0,®)  with  y.  It  is  this  relationship  which  is  the  basis  of  the  extreme 
value  limit  theorem  (see  Appendix  3),  and  in  that  sense  the  limit  theorem 
is  equivalent  to  the  asymptotic  relation  (11). 
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The  function  Dw(z)  and  its  differentiation  with  respect  to  w 


The  parabolic  cylinder  function  D  (z)  is  related  to  the  Gamma  func- 

w 

tion  and  confluent  hypergeometric  function  by  [12],  p.  324) 


CD 


1 

2 


where 


F(-l/2)  F  fl-w. 

r(-w/2)  1  ll  2  ' 


3 

2' 


(2) 


jF^v;  c;  z) 


Il£l  y  r,v*n)  z" 
r(v)  n=0  r(c+n)  n! 


and  this  representation  is  valid  for  all  complex  w  and  z.  The  factor 
2w^2exp (-z2/4)  is  entire  in  w  and  z  and  r'^v)  is  entire  in  v.  Hence, 
any  singularities  of  Dw(z)  in  the  complex  w  plane  must  be  associated  with 
jFjfv;  c;  z) .  We  note  that 

*  |v(v+l)  ...  (v*n-l)[  s  | v |  ( | v | + 1 )  ...  ( | v | +n-l) 


so  that  IjFjfv;  c;  z) |  s  |jFj(|v|;  c;  z)|.  Since  v(v+l) . . . (v*n-l)  is  entire 
in  v,  it  follows  from  the  inequality  above  and  the  Weierstrass  M-test  that 
jFj{v;  c;  z)  is  entire  in  v  and  z  and  Dw(z)  is  entire  in  w  and  z.  As 
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shown  in  [11],  D_g(-y)  has  zeros  only  for  s  real  and  negative. 

To  derive  equation  (2.10),  we  note  that 

(3)  f  e‘x2/4D  (-x)dx  -  /  e~x  /4D_g(x)dx 

-oo  -y  ** 

-7  -  fx[e'x2/4D.s.1(*)]dx  -  e"y2/4D-#-1(-y)  . 

The  second  equality  is  obtained  from  [12,  p.  327]  with  the  misprint  there 
corrected.  Equation  (2.10)  then  follows  from  (2.9),  and  the  representa¬ 
tion 


(4) 


D-S-l(-y)  7  en(y) 

'  n=l  5  *  Vy) 


arising  from  the  poles  at  -*n(y)  and  the  residues  there.  The  positivity  of 
Bn(y)  may  be  derived  from  the  complete  monotonicity  of  the  ergodic  exit  time 

E 

density  s^ft)  discussed  in  Section  2B.  Equation  (2.11)  is  obtained  from 
(2.8)  and  (2.10)  after  observing  from  (2.10)  that  a  (0)  *  1  implies  that 


(5) 


y,  '-*2'2  .  e'y2/2  7  «,<!-> 

L  /K  Jli  nil  Vy>  ■ 


To  verify  equation  (15)  and  thence  (16),  we  proceed  as  follows.  From  (2.8) 
and  (2.10)  we  have 


(6) 


Wx)>  •  /  - 


-x2/2 


i  *^2tT 


dx 


1  -y2/2  D-s-l('y) 

/Hi  *v.st-y) 
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It  follows  finally  that 


yfe(y)} 


i.e. , 


(10) 


yf,Mi 


2 


establishing  (2.15). 


Calculations  for  P [M (6 )  <  0]  *  Fg(0) 


The  p.d.f.  of  Tjq,  the  passage  time  from  x  <  0  to  x  =  0  is  known 
[11]  to  have  the  form 

2x(l*u)exp{-  i  x2/u) 

(1)  s(x,0,t)  »  - — — TTj - 


where  u  *  e2t  ••  1.  Hence  the  ergodic  exit  time  density  from  (-“,0)  is 


0  -x  /2 

(2)  s£(t)  *  2  /  —  s^c,0,T)dx 

-•  *2ri 


2  -1/2  2  -tn  -2t,  -1/2 

'  -  u  *  —  e  (1  -  e  ) 

7T  IT 


i.e.,  one  has  (2.19).  Eq  (2.20)  then  follows  and  the  coefficients  are  found 
from  the  identity 


(3) 


2 

"  F(j,  b;  b;  w2)  ,  any  b  , 


/l  -  w2 


2  w 


2n 


it  T (1/2)  F(1/2  +  n)  n! 


where  F  is  the  hypergeometric  function.  The  asymptotic  relation  (2.22)  is 
obtained  from  Stirling's  Formula. 
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Appendix  3 

Derivation  of  the  limit  theorem  for  extreme  values 

The  limit  theorem  for  M(6)  and  Ffi(y)  may  be  extracted  from  the  con¬ 
tents  of  Darling  and  Erdos  in  the  following  way. 

Def.  1.  Let  X(t)  be  the  stationary  0-U  process  as  described  in 
Section  1,  and  let  M(t)  be  the  maximum  process  as  described  there. 

Def.  2.  Let  T(a)  be  the  passage  time  of  X(t)  to  a,  i.e., 

T(a)  *  sup[t|x(t)  <  a,  0  s  t  s  t]. 

Def.  3.  Let  tj,  *  log  k. 

Def.  4.  Let  U  «=  max  X(t.)  *  M  (t  ) . 

-  n  ,  .  K  n 

isksn 

Def.  S.  Let  N(a),  1  s  k  s  n,  be  the  first  time  that  Un  reaches  the 
set  [a,®). 

Then  Darling  and  Erdos  show  that 


(a) 

P[Un  <  o]  *  P[N(o)  >  n]  «  P(log  N(a)  >  log  n] 

(L.3.2) 

(b) 

2. 

P[T(a)  >  {ea  Z/a)y]  ■*  e’y  ,  a  -*•  ® 

(L. 3.4) 

(c) 

p 

log  N(a)  -  2T(a)  •*  0  ,  a  -*•  ® 

(L.3.8) 

Hence, 

from  Def.  4,  and  (a) 

P[M  (tn)  <  o]  *  P[log  N(o)  >  log  n]  ,  all  a,  n  . 

•  «.  def  2y(a)y  .  .  .  ^2tT  a^/2 

Let  n  e  e  7  where  y(a)  *  -  e  , 


so  that 


P[M  (i  log  n  )  <  a]  «  P[log  N(a)  >  log  n  ]  .  all  a 

4  a  9/ 


From  (b) ,  (c) 


lim  P[M  (i  log  n  )  <  a}  *  lim  P[T(a)  >  y(a)y] 

..  4  a»;  _ 


e-y  for 


all  fixed  y  >  0 


If  we  now  set 


\  log  nQ  y  *  u(a)y  «=  0 


and  solve  (3)  for  a  as  a  function  of  y  and  0,  Eq.  (2)  becomes 


lim  P[M*(0)  <  a(y,0)]  -  e"y  ,  all  y  >  0  . 


1  2 
_1  j  a 

But  (3)  has  the  form  /lit  a-1  e  y  ■  0,  and  its  solution  for  a  is  given 
asymptotically  (see  below)  by 


.fy  6)  ~  c(0)  4  ^  Stf.?  -  I 
uy.oj  ciej  c(ej  2  c(0) 


,  0  -►  • 


where 


c(0)  «  /21og  0 
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and  this  is  the  form  of  the  limit  theorem  quoted  in  Section  3.  Darling 

* 

and  Erdos  state,  without  supplying  details,  that  M  (6)  and  M(6)  have  the 
same  asymptotic  behavior.  It  is  easier  to  verify  this  directly  from 
(3.5),  which  states  that 

(7)  P[M(0)  <  y]  =  P[M(0)  <  y]P[Ty  >  6]  . 

Hence 

(6)  P[M(u0yw)  <  y]  «  P[M(0)  <  y]P[(Ty/w0y)  >  w]  . 

But  E[Ty)  -  uQy  =  E[TQy]  as  y  -*•  »  and  P[T^/wy  >  w]  -►  e-w,  as  y  From 

(8)  we  then  have 

(9)  lim  P[M(6)  <  y(0 ,w) ]  -  e"w 

8-m» 

where  y(0,w)  is  the  solution  of 

e-y2/2w 

(10)  6  x  wv0y  ~  /2v  - - -  ,  y  -*•  - 

so  that  (5)  and  the  limit  theorem  follow. 

Asymptotic  solution 

The  equation  to  be  solved  may  be  written  as  v  x  ex/x  where  v  =  62/(2iry2), 
2 

x  *  a  .  Then  x  *  log  v  «•  log  x  ->  x  *  log  v  ♦  log  log  v  ♦  o(log  log  v) , 
v  -  and  a  »  /log  v  ♦  log  log  v  -  w  ♦  (log  w)/w  ♦  o((log  w)/w)  where 


w  =  /log  v  =  /21og6  -  log(2ity2)  =  c(6)  -  log (2ny2)/c (6) .  Hence 
(log  w)/w  *  (log  c(6))/c{0)  +  o (c  *),  and  finally  one  obtains  Eq.  (5), 
and  thence  Eq.  (4). 


r 
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1  TABLE  1 


c 

0.0 

0.0500 

0.1000 

0.1500 

0.2000 

0.2500 

THETA 

Y 

1.0000 

1.0013 

1.0050 

1.0113 

1.0202 

1.0317 

-9.00 

0.0000 

0,0000 

0.0000 

0.0000 

0.0000 

0.0000 

-3.90 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-3.80 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-3. TO 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-3.80 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-3.50 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

•3.90 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-3.30 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-3.20 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

•3.10 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-3.00 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

•2.90 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-2.80 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-2.70 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-2.60 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-2.50 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

•2,90 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-2.30 

0.0001 

0.0001 

0.0001 

0.0001 

0.0001 

0.0001 

-2.20 

0.0001 

0.0001 

0.0001 

0.0001 

0.0001 

0.0001 

-2.10 

0.0002 

0.0002 

0.0002 

0.0002 

0.0002 

0.0001 

-2.00 

0.0003 

0.0003 

0.0002 

0.0002 

0.0002 

0.0002 

-1.90 

0.0009 

0.0009 

0.0009 

0.0009 

0.0009 

0.0003 

-1,80 

0.0006 

0.0006 

0.0006 

0.0006 

0.0005 

0.0005 

-1,70 

0.0009 

0.0009 

0.0009 

0.0008 

0.0008 

0.0006 

-1,60 

0.0013 

0.0013 

0.0013 

0.0012 

0.0012 

0.0012 

-1.50 

0.0019 

0.0019 

0.0019 

0.0016 

o.ooie 

0.0017 

-1.90 

0.0027 

0.0027 

0.0027 

0.0026 

0.0026 

0.0025 

-1.30 

0.0038 

0.0036 

0.0038 

0.0037 

0.0036 

0.0035 

-1.20 

0.0059 

0.0053 

0.0053 

0.0052 

0.0051 

0.0050 

-1.10 

0.0079 

0.0079 

0.0073 

0.0072 

0.0071 

0.0069 

-1.00 

0.0101 

0.0101 

0,0100 

0.0099 

0.0097 

0.0095 

-0.90 

0.0137 

0.0137 

0.0136 

0.0139 

0.0132 

0.0129 

•0.80 

0.0183 

0.0182 

0,0181 

0.0179 

0.0176 

0.0172 

-0.70 

0.0291 

0.0290 

0,0239 

0.0236 

0.0233 

0.0228 

•0.60 

0.0319 

0.0313 

0,0311 

0.0306 

0.0309 

0.0298 

-0.50 

0.0909 

0.0903 

0,0901 

0.0397 

0.0392 

0.0385 

-0.90 

0.0513 

0.0513 

0.0510 

0.0505 

0.0999 

0.0991 

-0.30 

0.0696 

0.0699 

0.0691 

0.0636 

0.0629 

0.0619 

-0.20 

0.0602 

0.0601 

0,0798 

0.0791 

0.0763 

0.0772 

-0.10 

0.098b 

0.0985 

0,0961 

0.0979 

0.0969 

0.0952 

2 


TABLE  I 


C 

THETA 

V 

4.00 

3.90 

3.80 

3.70 

3.60 

3.50 

3.40 

3.30 

3.20 
3. tO 
3.00 

2.90 

2.80 

2.70 

2.60 

2.50 

2.40 

2.30 

2.20 
2.10 
2.00 

1.90 
1.80 

1.70 
1.60 

1.50 

1.40 

1.30 
1.20 
1.10 
1.00 
0.90 
0.80 
0.70 
0.60 
0.50 
0.40 
0.30 
0.20 
0.10 


0.3000 

1.0460 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0001 

0.0001 

0.0001 

0.0002 

0.0003 

0.0005 

0.0008 

0.0011 

0.0017 

0.0024 

0.0034 

0.0048 

0.0067 

0.0092 

0.0125 

0,0168 

0.0223 

0.0291 

0.0377 

0.0481 

0.0608 

0.0759 

0.0936 


0.3500 

1.0632 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0001 

0.0001 

0.0001 

0.0002 

0.0003 

0.0005 

0.0007 

0.0011 

0.0016 

0.0023 

0.0033 

0.*>046 

0.0065 

0.0089 

0.0121 

0.0163 

0.0216 

0.0263 

0.0367 

0.0470 

0.0594 

0.0743 

0.0918 


0.4000 

1.0833 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0001 

0.0001 

0.0002 

0.0003 

0.0004 

0.0007 

o.ooto 

C.0015 

0.0022 

0.003! 

0.0044 

0.0062 

0.0086 

0.0117 

0.0157 

0.0209 

0.0274 

0.0356 

0,0457 

0.0579 

0.0725 

0.0698 


0.4500 

1.1066 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0001 

0.0001 

0,0002 

0.0003 

0.0004 
0,0006 
0.0010 
0.0014 
0.0021 
0.0030 
0.0042 
0.0059 
0.0082 
0.011 1 
0.0151 
0.0201 
0.0264 
0.0344 
0.0442 
0.056? 
0.0705 
0.0675 


0.5000 

1.1331 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0001 

0.0001 

0.0002 

0.0002 

0.0004 

0.0006 

0.0009 

0.0013 

0.0019 

0.0028 

0.0039 

0.0056 

0.0077 

0.0106 

0.0143 

0.0192 

0,0253 

0.0331 

0.0426 

0.0543 

0.0683 

0.0849 


0.5500 
1 .1633 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

o.ooco 

0.0000 

0.0000 

0.0001 

0.0001 

0.0001 

0.0002 

0.0003 

0.0005 

0,0008 

0.0012 

0.0018 

0.0026 

0.0037 

0.0052 

0.0073 

0.0100 

0.0136 

0.0182 

0.024? 

0.0316 

0.0409 

0.0522 

0.0658 

0.0821 
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3  TABLE  1 


c 

0.6000 

0.6500 

0.7000 

0.7500 

0.6000 

0.6500 

T«ETA 

Y 

1.1972 

1.2352 

1.2776 

1.3296 

1.3771 

1.9351 

-a. oo 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-3.90 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-3.60 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-3. TO 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-3.60 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-3. SO 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-3.90 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-3.30 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-3.20 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-3.10 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-3.00 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-2.90 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-2.80 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-2.70 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-2.60 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-2.50 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-2.9  0 

0.0000 

A. 0000 

0.0000 

0.0000 

0.0000 

0.0000 

-2.30 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-2.20 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-2.10 

0.0001 

0.0001 

0.0001 

0.0000 

0.0000 

0.0000 

-2.00 

0.0001 

0.0001 

0.0001 

0.0001 

0.0001 

0.0001 

-1.90 

0.0002 

0.0002 

0.0002 

0.0001 

0.0001 

0.0001 

-1.80 

0.0003 

0.0003 

0.0002 

0.0002 

0.0002 

0.0001 

-1.70 

0.0005 

0.0009 

0.0009 

0.0003 

0.0003 

0.0002 

-1.60 

0.0007 

0.0007 

0.0006 

0.0005 

0.0009 

0.0009 

-1.50 

0.0011 

0.0010 

0.0009 

0.0006 

0.0007 

0.0006 

-1.90 

0.0016 

0.0015 

0.0013 

0.0012 

0.0010 

0.0009 

-1.30 

0.0029 

0.0022 

0.0020 

0.0017 

0.0015 

0.0013 

-1.20 

0.0039 

0.0031 

0.0026 

0.0025 

0.0023 

0.0020 

-1.10 

0.0096 

0.0095 

0.0091 

0.0037 

0.0033 

0.0029 

-1.00 

0.0068 

0.0063 

0.0058 

0.0052 

0.0097 

0.0092 

•0.90 

0.0099 

0.0087 

0.0060 

0.0073 

0.0066 

0.0059 

-0.60 

0.0128 

0.0119 

0.0111 

0.0102 

0.0092 

0.0083 

-0.70 

0.0172 

0.0161 

0.0150 

0.0139 

0.0127 

0.0115 

-0.60 

0.0229 

0.0215 

0.0201 

0.0187 

0.0172 

0.0157 

-0.50 

0.0301 

0.0269 

0.0267 

0.0299 

0.0230 

0.0211 

-0.90 

0.0390 

0.0370 

0.0399 

0.0327 

0.0309 

0.0261 

-0.30 

0.0999 

0.0975 

0.0950 

0.0929 

0.0396 

0.0366 

•0.20 

0.0632 

0.0609 

0.0579 

0.0592 

0.0510 

0.0976 

-0.10 

0.0790 

0.0758 

0.0723 

0.0666 

0.0697 

0.0607 

9 
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■ 

« 

;  ' 

c 

0,9000 

0.9500 

;  • 

THETA 

Y 

1.0903 

1.5703 

1 

•9.00 

0,0000 

0.0000 

, 

-5,90 

0.0000 

0.0000 

-3.80 

0.0000 

0.0000 

-3.70 

0.0000 

0.0000 

-3.60 

0.0000 

0.0000 

i 

-3.50 

0.0000 

0.0000 

-3.00 

0.0000 

0.0000 

-3.30 

0.0000 

0.0000 

-3.20 

0.0000 

0.0000 

-3.10 

0.0000 

0.0000 

-3.00 

0.0000 

0,0000 

-2.90 

0.0000 

0,0000 

-2.80 

0.0000 

0,0000 

{ 

-2.70 

0.0000 

0,0000 

-2.60 

0.0000 

0.0000 

\ 

-2.50 

0.0000 

o.oooc 

\ 

-2.«0 

0.0000 

0.0000 

-2.30 

0.0000 

0.0000 

-2.20 

0.0000 

0,0000 

-2.10 

0.0000 

0.0000 

-2.00 

0.0000 

0.0000 

-1.90 

0.0001 

0.0001 

-1.80 

0.0001 

0.0001 

-1.70 

0.0002 

0.0002 

-1.60 

0.0003 

0.0003 

-1.50 

0.0005 

0.0009 

| 

-1.00 

0.0008 

0.0006 

-1.30 

o.oon 

0.0010 

1 

-1.20 

0.0017 

0.0015 

1 

-1.10 

0.0025 

0.002? 

i 

-1.00 

0.0037 

0.0032 

-0.90 

0.0053 

0,0096 

I 

-0.80 

0.0079 

0.0065 

t 

-0.70 

0.0103 

0.0092 

i 

•0.60 

0.0192 

0.0127 

! 

-0.50 

0.0192 

0.0173 

-0.00 

0.0257 

0.0233 

-0.30 

0.0339 

0,0310 

-0.20 

0.0991 

0.0905 

-0.10 

0.0566 

0.0523 

TA8LE  1 


1.0000 

1.0500 

1.1000 

1.1500 

1.6987 

1.7359 

1.6313 

1.9372 

0.0000 

0,0000 

0.0000 

0.0000 

0.0000 

0,0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

o.oooo 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0001 

0.0001 

0.0000 

0.0000 

0.0001 

0.0001 

0.0001 

0.0001 

0.0002 

0.0002 

0.0001 

0.0001 

0.0003 

0.0003 

0.0002 

0.0002 

0.0005 

0.0009 

0.0003 

0.0002 

0.0008 

0.0007 

0.0005 

0.0009 

0.0012 

0.0010 

0.0006 

0.0006 

o.ooie 

0.0015 

0.0012 

0.0010 

0.0027 

0.0023 

0.0019 

0.0015 

0,0090 

0.0039 

0.0028 

0.0023 

0,0057 

0.0099 

0.0091 

0.0039 

0.0061 

0.0070 

0.0059 

0.0050 

0,0112 

0.0098 

0.0085 

0.0072 

0.0155 

0.0136 

0.0118 

0.0102 

0.0210 

0.0166 

0.0169 

0.0192 

0.0260 

0.0251 

0.0222 

0.0199 

0.0369 

0.0333 

0.0298 

0.0263 

0.0960 

0.0937 

0.0393 

0.0351 
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5 


C 

TmET* 

Y 

a. 00 

3. 90 

3. so 

3.70 

3.60 

3. 50 

3.40 

3.30 

3.20 

3.10 

3.00 

2.90 

2. 50 

2.70 

2.60 

2. 50 

2.40 

2.30 

2.20 

2.10 
2.00 

1.90 
1.80 

1.70 
1.60 

1.50 

1.40 

1.30 
1.20 
1.10 
1.00 
0.90 
0.80 
0.70 
0.60 
O.SO 
0.40 
0.30 
0.20 
0.10 


TABLE  1 

1.2000 

1.2500 

1.3000 

2.0S44 

2.1842 

2.3280 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

O.V/OOO 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0001 

0.0000 

0.0000 

0.0001 

0.0001 

0.0001 

0.0002 

0.0001 

0.0001 

0.0003 

0.0002 

0.0002 

o.ooos 

0.0004 

0.0003 

0.0008 

0.0006 

0.0004 

0.0012 

0.0009 

0.0007 

0.0018 

0.0014 

0.0011 

0.0028 

0.0022 

0.0017 

0.0041 

0.0033 

0.0026 

0.0060 

0.0049 

0.0039 

0.0086 

0.0071 

0. 00S6 

0.0121 

0.0101 

0.0083 

0.0168 

0.0142 

0.0119 

0.0229 

0.0197 

0.0166 

0,0309 

0.0268 

0.0229 

1.3500 

i.4000 

1.4500 

2.4674 

2. 6645 

2.8612 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

o.oooo 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0001 

0.0000 

0.0000 

0.0001 

0.0001 

0.0000 

0.0002 

0.0001 

0.0001 

0.0003 

0.0002 

0.0001 

0.0005 

0.0004 

0.0002 

0.0006 

0.0006 

0.0004 

0.0013 

0.0009 

0.0007 

0.0020 

0,0015 

0.0011 

0.0030 

0.0023 

0.0017 

0.0046 

0.0035 

0.0027 

0.0067 

0.0053 

0.0041 

0.0097 

0.0078 

0.0061 

0.0136 

0.0112 

0.0069 

0.0193 

0,0159 

0.0129 
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c 

1.5000 

1.5500 

THETA 

V 

3.0802 

3.3243 

-4.00 

0.0000 

0.0000 

-1.90 

0.0000 

0.0000 

-3.80 

0.0000 

0.0000 

-1.70 

0.0000 

0.0000 

•1.60 

0.0000 

0.0000 

-1.50 

0.0000 

0,0000 

-3.40 

0.0000 

0.0000 

-3.10 

0.0000 

0.0000 

-3.20 

0.0000 

0.0000 

-3.10 

0.0000 

0.0000 

-3.00 

0.0000 

0.0000 

-2.90 

0.0000 

0.0000 

-2.80 

0.0000 

0.0000 

-2.70 

0.0000 

0.0000 

•2.60 

0.0000 

0.0000 

-2.50 

0.0000 

0.0000 

-2.40 

0.0000 

0.0000 

-2.30 

0.0000 

0.0000 

-2.20 

0.0000 

0.0000 

-2.10 

0.0000 

0.0000 

-2.00 

0.0000 

0.0000 

-1.90 

0.0000 

0,0000 

-1.80 

0.0000 

0.0000 

-1.70 

0.0000 

0.0000 

-1.60 

0.0000 

0.0000 

•1  .50 

0.0000 

0.0000 

-1.40 

0.0000 

0.0000 

-1.30 

0.0000 

0.0000 

-1.20 

0.0000 

0,0000 

-I. to 

0.0001 

0.0001 

-1.00 

0.0002 

0,0001 

-0.90 

0.0003 

0.0002 

-0.80 

0.0005 

0.0003 

-0.70 

0.0007 

0.0005 

-0.60 

0.0012 

0.0008 

-0.50 

0.0019 

0.0014 

•0.40 

0.0030 

0.0022 

•0.30 

0,0046 

0.0034 

•0.20 

0.0069 

0.0052 

-0.10 

0.0101 

0.0078 

table  1 

1.6000 

1 .6500 

1.7000 

3.5966 

3.9011 

4.2419 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0,0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0,0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

o.cooo 

0.0000 

0,0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0,0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0001 

0.0000 

0.0000 

0,0001 

0.0001 

0.0000 

0.0002 

0.0001 

0.0001 

0.0003 

0.0002 

0.0001 

0.0005 

0.0003 

0.0002 

0,0009 

0.0006 

0.0004 

0.0015 

0.0010 

0.0006 

0.0024 

0.0016 

0.0011 

0.0038 

0.0026 

0.0018 

0.0058 

0.0042 

0.0029 

1 .7500 
4.62«0 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0001 

0.0001 

0.0002 

O.OOOU 

0.0007 

0.0011 

0.0019 


-40- 


7 


c 

1.8000 

1.8500 

THETA 

Y 

5.0531 

5.5359 

-6.  00 

0.0000 

0.0000 

-3.90 

0.0000 

0.0000 

-3.80 

0.0000 

0.0000 

-3.70 

0.0000 

0.0000 

-3.60 

0.0000 

0.0000 

-3.50 

0.0000 

0.0000 

-3.60 

0.0000 

0.0000 

-3.30 

0.0000 

0.0000 

-3.20 

0.0000 

0.0000 

-3.10 

0.0000 

0.0000 

-3.00 

0.0000 

0.0000 

-2.90 

0.0000 

0.0000 

-2.80 

0.0000 

0.0000 

-2.70 

0.0000 

0.0000 

-2.60 

0.0000 

0.0000 

-2.50 

0.0000 

0.0000 

-2.60 

0.0000 

0.0000 

-2.30 

0.0000 

0.0000 

-2.20 

0.0000 

0.0000 

-2.10 

0.0000 

O.OOQO 

-2.00 

0.0000 

0.0000 

-1.90 

0.0000 

0.0000 

-1.80 

0.0000 

0.0000 

-1.70 

0.0000 

0.0000 

-1.60 

0.0000 

0.0000 

-1.50 

0.0000 

0.0000 

-l.«0 

0.0000 

0.0000 

-1.30 

0.0000 

0.0000 

-1.20 

0.0000 

0.0000 

•1.10 

0.0000 

0.0000 

-1.00 

0.0000 

0.0000 

-0.90 

0.0000 

0.0000 

-0.80 

0.0000 

0.0000 

-0.70 

0.0000 

0.0000 

-0.60 

0.0001 

0.0000 

-0.50 

0.0001 

0.0001 

•0.60 

0.0002 

0.0001 

-0.30 

o.oooa 

0.0002 

•0.20 

0.0007 

0.0006 

-0.10 

0.0012 

0,0007 

TABLE  1 

1.9000 

1.9500 

2.0000 

6.0600 

6.6963 

7.3691 

0.0000 

0.0 

0.0 

0.0000 

0,0000 

0.0 

0.0000 

0.0000 

0.0 

0.0000 

0.0000 

0.0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0,0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0001 

0.0000 

0.0000 

0.0001 

0.0000 

0.0000 

0.0002 

0.0001 

0.0000 

0.0006 

0.0002 

0.0001 

2.050  0 
6. 1766 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0,0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


r 

r- 

-41- 
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TABLE  1 

c 

2.1000 

2.1500 

2.2000 

2.2500 

2.3000 

2.3500 

theta 

y 

9.0703 

10.0870 

11.2959 

12,5692 

19.0839 

15.8196 

}j 

-a. oo 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-3.90 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-3.80 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

, 

-3.70 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

y 

-3.60 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-3.50 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-3.90 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

?  ■ 

-3.30 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-3.20 

0,0 

0.0 

0.0 

0.0 

0.0 

0.0 

-3.10 

0.0000 

0.0 

0.0 

0.0 

0.0 

0.0 

-3.00 

0.0000 

0.0 

0.0 

0.0 

0.0 

0.0 

-2.90 

0.0000 

0.0 

0.0 

0,0 

0.0 

0.0 

-2.80 

0.0000 

0.0000 

0.0 

0.0 

0.0 

0.0 

-2.70 

0.0000 

0.0000 

0.0 

0.0 

0.0 

0.0 

-2.60 

0.0000 

0.0000 

0.0 

0.0 

0.0 

0.0 

-2.50 

0.0000 

0.0000 

0.0000 

0.0 

0.0 

0.0 

-2.90 

0.0000 

0.0000 

0.0000 

0.0 

o,  c 

0.0 

-2.30 

0.0000 

0.0000 

0.0000 

0.0000 

0.0 

0.0 

. 

-2.20 

0.0000 

0.0000 

0.0000 

0,0000 

0.0 

0.0 

-2.10 

0.0000 

0.0000 

0.0000 

0.0000 

0.0 

0.0 

•2.00 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0 

-1.90 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0 

-1.80 

0.0000 

0.0000 

0.0000 

0,0000 

0.0000 

0.0000 

-1.70 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-1.60 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-1.50 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

. 

‘ 

-1.90 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-1.30 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-1.20 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

•1.10 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-1.00 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.90 

0.0000 

0.0000 

0.0000 

0,0000 

0.0000 

0.0000 

-0.80 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.70 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0,60 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.50 

0.0000 

0.0000 

0.0000 

0,0000 

0.0000 

0.0000 

1 

-0.90 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

: 

-0.30 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.20 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

L 

•  -0.10 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

: 

; 

j 
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9  TABLE  l 


c 

0.0 

0.0500 

0.1000 

0.1500 

0.2000 

0.2500 

THETA 

Y 

1.0000 

1.0013 

1.0050 

1.0113 

1.0202 

1.0317 

0.0 

0.1144 

0.1198 

0.1193 

0.1185 

0.1174 

0.1160 

0.10 

0.1442 

0.1490 

0.1435 

0.1426 

0.1414 

0.1398 

0.20 

0.1716 

0.1714 

0.1708 

0.1698 

0.1685 

0.1668 

0.30 

0.2020 

0.2018 

0.2012 

0.2001 

0.1967 

0.1968 

0.40 

0.2354 

0.2352 

0.2345 

0.2334 

0.2319 

0.2299 

O.SO 

0.2716 

0.2714 

0.2707 

0.2695 

0.2679 

0.2656 

0.60 

0.3104 

0.3101 

0.3094 

0.3082 

0, 3065 

0.3043 

0.70 

0.3513 

0.3510 

0.3503 

0.3490 

0.3472 

0.3450 

o.eo 

0.3434 

0.3936 

0.3929 

0.3916 

0.3698 

0.3875 

0.40 

0.4378 

0.4375 

0.4368 

0.4355 

0.4337 

0.4314 

1.00 

0.4824 

0.4821 

0.4814 

0.4801 

0.4784 

0.4761 

1.10 

0.5271 

0.5269 

0.5262 

0.5249 

0.5232 

0.5210 

1.20 

0.5715 

0.5713 

0.5706 

0.5694 

0.5677 

0.5656 

1.30 

0.6144 

0.6147 

0.6140 

0.6129 

0.6113 

0.6093 

1.40 

0.6564 

0.6567 

0.6560 

0.6550 

0.6535 

0.6516 

1.50 

0.6464 

0.6967 

0.6961 

0.6952 

0.6938 

0.6920 

1.60 

0.7347 

0.7345 

0.7340 

0.7331 

0.7318 

0.7302 

1.70 

0.7648 

0.7697 

0.7692 

0.7684 

0.7672 

0.7658 

l.BO 

0.8022 

0.8020 

0.8016 

0.8009 

0.7999 

0.7965 

1.40 

0.8316 

0.8315 

0.8311 

0.8304 

0.8295 

0.6284 

2.00 

0.8580 

0.8579 

0.8575 

0.8570 

0.6562 

0.8552 

2.10 

0.8814 

0.8813 

0.8810 

0.6805 

0.8796 

0.8789 

2.20 

0.4014 

0.9018 

0.9015 

0.9011 

0.9006 

0.8998 

2.30 

0.4146 

0.9196 

0.9194 

0.9190 

0.9165 

0.9179 

2.40 

0.4148 

0.9348 

0.9346 

0.9343 

0.9339 

0.9334 

2.50 

0.4477 

0.9476 

0.9475 

0.9473 

0.9469 

0.9465 

2.60 

0.4584 

0.9584 

0.9563 

0.9581 

0.9578 

0.9574 

2.70 

0.4673 

0.9672 

0.9672 

0.9670 

0.9666 

0.9665 

2. 80 

0.4745 

0.9745 

0.9744 

0.9743 

0.9741 

0.9739 

2.40 

0.4803 

0.9803 

0,9803 

0.9802 

0.9800 

0.9799 

3.00 

0.4850 

0.9850 

0.9849 

0.9649 

0.9848 

0.9646 

3.10 

0.4887 

0.9887 

0.9866 

0.9886 

0.9865 

0,9684 

3.20 

0.4415 

0.9915 

0.9915 

0.9914 

0.9914 

0.9913 

3.30 

0.4437 

0.9937 

0.9937 

0.9937 

0.9936 

0.9936 

3.40 

0.4454 

0.9954 

0.9954 

0.9954 

0.9953 

0.9953 

3.50 

0.4467 

0.9967 

0.9967 

0.9966 

0.9966 

0.9966 

3.60 

0.4476 

0.9976 

0.9976 

0.9976 

0.9976 

0.9975 

3.70 

0.4483 

0.9983 

0.9983 

0.9963 

0.9983 

0.9983 

3.80 

0.4488 

0.9988 

0.9988 

0.9988 

0.9986 

0.9968 

3.40 

0.4442 

0.9992 

0.9992 

0.9992 

0.9992 

0.9991 

4.00 

0.4494 

0.9994 

0.9994 

0.9994 

0.9994 

0.9994 

4.10 

0.4446 

0.9996 

0.9996 

0.9996 

0.9996 

0.9996 

4.20 

0.4997 

0.9997 

0.9997 

0.9997 

0.9997 

0.9997 

4.30 

0.4998 

0.9998 

0.9998 

0.9998 

0.9998 

0,9998 

4,40 

0.9999 

0,9999 

0.9999 

0.9999 

0.9999 

0.9999 

4.50 

0.9494 

0.9999 

0.9999 

0,9999 

0.9999 

0.9999 
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10  table  1 


c 

0.3000 

0.3500 

0.0000 

0.0500 

0.5000 

THETA 

r 

1.0060 

1.0632 

1.0833 

1.1066 

1.1331 

0.0 

0. 1 103 

0.1122 

0.1099 

0.1073 

0.10OO 

0.10 

0.1379 

0.1357 

0.1331 

0.1301 

0.1266 

0.20 

0.1607 

0.1622 

0.1593 

0.1561 

0. 152o 

0.30 

0.1905 

0.1918 

0.1867 

0.1652 

0.1812 

0.00 

0.227O 

0.2205 

0.2212 

0.2170 

0.2131 

0.50 

0.2632 

0.2601 

0.2566 

0.2525 

0.2060 

0.60 

0.3016 

0.2980 

0.2907 

0,2900 

0.2857 

0.70 

0.3022 

0.3369 

0.3351 

0.3307 

0.3258 

0.80 

0.38O7 

0.3610 

0.3775 

0.3730 

0.3681 

0.90 

0.9286 

0.O252 

0.O213 

0.0169 

0.0118 

1.00 

0.0733 

0.0700 

0 .066 1 

0.0617 

0.0567 

1.10 

0.5183 

0.5150  . 

0.5112 

0.5069 

0.5020 

1.20 

0.5630 

0.5598 

0.5562 

0.5520 

0.5072 

1.30 

0.6068 

0.6038 

0.6003 

0.5963 

0.5918 

1.00 

0.6092 

0.6A60 

0 .6031 

0  «  6390 

0.6351 

1.50 

0.6696 

0.6872 

0.6602 

0.6806 

0.6767 

1.60 

0.7282 

0.7258 

0.7230 

0.7197 

0.7161 

1.70 

0.7600 

0.7618 

0.7592 

0.7563 

0.7530 

1.60 

0.7969 

0.7950 

0.7927 

0.7901 

0.7871 

1.90 

0.8269 

0.8252 

0.6232 

0.6209 

0.8182 

2.00 

0.6539 

0.6520 

0.6506 

0.6086 

0.8063 

2.10 

0.8779 

0.8766 

0.8751 

0.8733 

0.8713 

2.20 

0.8989 

0.8970 

0.8965 

0.8950 

0.8933 

2.30 

0.9171 

0.9162 

0.9151 

0.9136 

0.9120 

2.00 

0.9327 

0.9320 

0.9311 

0.9300 

0.9268 

2.50 

0 .9060 

0 .9053 

0.9006 

0.9O37 

0.9027 

2.60 

0.9570 

0.9565 

0.9559 

0.9552 

0.9500 

2.70 

0.9662 

0.9657 

0,9653 

0.9607 

0.9600 

2.80 

0.9736 

0.9733 

0.9729 

0.9725 

0.9719 

2.90 

0.9796 

0 .9790 

0.9791 

0.9787 

0.9763 

3.00 

0.9605 

0.9603 

0,9800 

0.9636 

0.9830 

3.10 

0.9883 

0.9861 

0,9079 

0.9677 

0.9875 

3.20 

0.9912 

0.9911 

0.9910 

0.9908 

0.9906 

3.30 

0.9935 

0.9930 

0.9933 

0.9932 

0.9930 

3.00 

0.9952 

0.9952 

0.9951 

0.9950 

0.9909 

3.50 

0.9965 

0.9965 

0.9960 

0.9960 

0.9963 

3.60 

0.9975 

0.9975 

0.9970 

0.9970 

0.9973 

3.70 

0.9982 

0.9982 

0.9962 

0.9981 

0.9981 

3.80 

0.9988 

0.9987 

0.9987 

0.9987 

0.9967 

3.90 

0.9991 

0,9991 

0.9991 

0.9991 

0.9991 

9.00 

0.9990 

0.9990 

0,9990 

0.9990 

0 .9990 

9.10 

0.9996 

0,9996 

0.9996 

0.9996 

0.9996 

9.20 

0.9997 

0.9997 

0.9997 

0.9997 

0.9997 

9.30 

0.9998 

0.9998 

0.9996 

0.9998 

0.9998 

9.90 

0.9999 

0.9999 

0,9999 

0.9999 

0.9999 

9.50 

0.9999 

0,9999 

0.9999 

0.9999 

0.9999 

o.ssoo 

1.1633 


0.1012 
0.1232 
o. ioeo 
0.1760 
0.2000 
0.2«30 

0.280O 
0.3200 
0.3625 
0.0062 
0.0511 
0.0065 
0.5019 
0.5867 
0.6303 
0.6722 
0.7120 
0.7O93 
0.7638 
0.0153 
0.8037 
0.8601 
0.0919 
0.9108 
0.9275 
0.9916 
0.9535 
0.9633 
0,9710 
0.9779 
0.9831 
0.9872 
0.9900 
0,9929 
0 . 9990 
0.9962 
0.9973 
0.9981 
0.9986 
0.9991 

0.9993 

0.9996 

0.9997 

0.9998 

0.9999 

0.9999 
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c 

0.6000 

0.6500 

0.7000 

0.7500 

0.0000 

0.0500 

theta 

Y 

1.1072 

1 .2352 

1,2776 

1.3240 

1.3771 

1.4351 

0.0 

0.0077 

0,0039 

0.0099 

0.0057 

0.0012 

0.0765 

0.10 

0.1103 

0.1151 

0.1105 

0.1057 

0.1006 

0.0952 

0.20 

0.1441 

0.1304 

0.1343 

0.1209 

0.1231 

0.1171 

o.so 

0.1721 

0.1660 

0.1613 

0.1553 

0.1489 

0.1422 

0.  AO 

0.2032 

0.1076 

0.1915 

0.1050 

0, 1701 

0.1707 

o.so 

0.2375 

0.2315 

0.2250 

0.2100 

0.2105 

0.2025 

0.60 

0.2747 

0.2603 

0.2615 

0.2540 

0.2461 

0.2376 

0.70 

0.3144 

0.3070 

0.3007 

0.2929 

0.2846 

0.2757 

o.eo 

0.3563 

0.3406 

0.3422 

0.3343 

0.3257 

0.3164 

0.00 

0.4000 

0.3032 

0.3057 

0.3777 

0.3609 

0.359U 

1.00 

0.4440 

0.4301 

0.4307 

0.4225 

0.4137 

0,4042 

1.10 

0.4004 

0.4037 

0.4764 

0.4604 

0.4507 

0.4502 

1.20 

0.5360 

0.5205 

0.5224 

0.5145 

0.5060 

0.4968 

1.30 

0.5010 

0.5740 

0.5679 

0.5604 

0.5522 

0.5433 

1.40 

0.6250 

0.6100 

0.6125 

0.6054 

0.5976 

0,5891 

1.50 

0.6672 

0.6617 

0.6556 

0.6409 

0.6416 

0,6336 

1.60 

0.7074 

0.7023 

0.6967 

0.6905 

0.6037 

0.6763 

1.70 

0.7451 

0.7405 

0.7353 

0.7297 

0.7235 

0.71b7 

i.eo 

0.7000 

0.7750 

0.7712 

0.7662 

0,7606 

0,7544 

1.00 

0.6110 

0.0003 

0.0042 

0.7996 

0.7947 

0.7892 

2.00 

0.0400 

0.0376 

0.0340 

0.0300 

0.0256 

0.6208 

2.10 

0.0666 

0.0637 

0.0606 

0.0572 

0.0534 

0.8492 

2.20 

0.0002 

0.0060 

0.0042 

0.0012 

0.0700 

0.0744 

2.30 

0.0000 

0 . 0070 

0.9047 

0.0022 

0.0995 

0.8965 

2.40 

0.0260 

0.0243 

0.9224 

0.0203 

0.0180 

0.9155 

2.50 

0.0404 

0.0390 

0.9374 

0.0357 

0.0338 

0.9317 

2.60 

0.0525 

0.0513 

0.0501 

0.0407 

0.0471 

0.9454 

2.70 

0.0625 

0.0616 

0.0606 

0.0505 

0.0562 

0.9568 

2.00 

0.0707 

0.0700 

0.0692 

0.9603 

0.0673 

0.9662 

2.00 

0.0774 

0.0760 

0.9762 

0.0755 

0.0747 

0.9738 

3.00 

0.0027 

0.0023 

0.0010 

0.0812 

0.0606 

0.9799 

3.10 

0.0060 

0.0066 

0.9062 

0,0050 

0.0853 

0.9648 

3.20 

0.0002 

0.0099 

0.9096 

0,0803 

0.0090 

0.9886 

3.30 

0.0027 

0.9925 

0.0923 

0.0921 

0.0918 

0.9915 

3.40 

0.0047 

0.9945 

0.0944 

0.9942 

0.9940 

0.9938 

3.50 

0.0061 

0.9960 

0.9959 

0.9950 

0.9956 

0.9955 

3.60 

0.0072 

0.9971 

0.9971 

0.9970 

0.9969 

0.9967 

3.70 

0.0900 

0.9900 

0.9979 

0.9970 

0.0970 

0.9977 

3.00 

0.0006 

0.9906 

0.9905 

0.9905 

0.9964 

0.9984 

3.00 

0.0000 

0.0990 

0.9990 

0,0909 

0.9989 

0.9909 

4.00 

0.0003 

0.9993 

0.0993 

0.9993 

0.9992 

0.9992 

4.10 

0.0005 

0.9995 

0,0095 

0.9995 

0.0995 

0.9995 

4.20 

0.0007 

0.0997 

0,9997 

0.9997 

0.9997 

0.9996 

4.30 

0.0000 

0.0990 

0.0090 

0.9990 

0.0990 

0.9998 

4.40 

0.0000 

0.0999 

0,9009 

0.9999 

0.9990 

0.9998 

4.50 

0.0000 

0.0999 

0.9099 

0.9999 

0,0999 

0.9999 

Ip- - “ 
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4 

U 

c 

0.9000 

0.9500 

table  i 

1.0000 

1.0500 

1.1000 

1.1500 

THETA 

Y 

1.4993 

1.5703 

1.6487 

1.7354 

1.8313 

1.9372 

j 

0.0 

0.0717 

0.0667 

0.0616 

0.0564 

0.0512 

0.0460 

\ 

0.10 

0.0897 

0.0839 

0.0780 

0.0719 

0.0656 

0.0596 

0.20 

0.1106 

0.1042 

0.0974 

0.0904 

0.0832 

0.0760 

0.30 

0.1351 

0.1277 

0.1201 

0.1121 

0.1040 

0.0957 

0.40 

0.1629 

0.1547 

0.1462 

0.1373 

0.1262 

0.1166 

0.50 

0.1941 

0.1852 

0.1758 

0.1660 

0,1559 

0.1454 

0.60 

0.2265 

0.2190 

0.2089 

0.1983 

0.1673 

0.1756 

0.70 

0.2661 

0.2560 

0.2453 

0.2340 

0.2222 

0.2096 

o.eo 

0.3065 

0.2960 

0.2646 

0.2729 

0.2604 

0.2472 

o.oo 

0.3493 

0.3365 

0.3269 

0.3146 

0.3016 

0.2679 

1.00 

0.3940 

0.3630 

0.3713 

0.3566 

0.3454 

0.3313 

1.10 

0.4400 

0.4291 

0.4173 

0.4047 

0.3913 

0.3770 

1.20 

0.4868 

0.4760 

0.4644 

0.4519 

0.4366 

0.4243 

1.30 

0.5336 

0.5231 

0.5118 

0.4997 

0.4666 

0.4726 

l.«0 

0.5799 

0.5696 

0.5590 

0.5473 

0.5347 

0.5211 

1.50 

0.6249 

0.6155 

0.6052 

0.5941 

0.5621 

0.5692 

1.60 

0.6662 

0.6594 

0.6499 

0.6395 

0.6263 

0.6161 

1.70 

0.7093 

0.7012 

0.6925 

0.6829 

0.6725 

0.6612 

1.60 

0.7478 

0.7404 

0.7325 

0.7236 

0.7143 

0,7040 

1.90 

0.7832 

0.7767 

0.7696 

0.7618 

0.7533 

0.7440 

2.00 

0.6156 

0.6096 

0.6035 

0.7966 

0.7891 

0.7808 

2.10 

0.8447 

0.6396 

0.6341 

0.6261 

0.6215 

0.8143 

2.20 

0.8705 

0.8662 

0.8614 

0.6563 

0.6506 

0-.8443 

2.30 

0.6931 

0.6695 

0.8654 

0.6610 

0.6762 

0.8709 

2.«0 

0.9127 

0.9096 

0.9063 

0.9026 

0.8965 

0.8940 

2. SO 

0.9294 

0.9269 

0.9241 

0.9210 

0.9177 

0.9140 

2.60 

0.9436 

0.9415 

0.9392 

0.9367 

0.9339 

0.9309 

j 

2.70 

0.9553 

0.9536 

0.9518 

0.9496 

0.9475 

0.9451 

| 

2.60 

0.9650 

0.9637 

0.9622 

0.9606 

0.9566 

0.9568 

2.90 

0.9729 

0.9718 

0.9707 

0.9694 

0.9660 

0.9664 

3.00 

0.9792 

0,9784 

0.9775 

0.9765 

0.9754 

0.9742 

3.10 

0.9842 

0.9836 

0.9829 

0.9621 

0.9813 

0.9804 

3*20 

0.9861 

0.9877 

0.9871 

0.9666 

0.9859 

0.9852 

3.30 

0.9912 

0.9906 

0.9904 

0.9900 

0.9895 

0.9890 

3.40 

0.9935 

0.9933 

0.9930 

0.9926 

0.9923 

0.9919 

3.50 

0.9953 

0.9951 

0.9949 

0.9946 

0.9944 

0.9941 

3.60 

0.9966 

0.9965 

0.9963 

0.9961 

0.9959 

0.9957 

3.70 

0.9976 

0.9975 

0.9974 

0.9972 

0.9971 

0.9970 

3.60 

0.9983 

0.9982 

0,9981 

0.9981 

0.9980 

0.9979 

3.90 

0.9988 

0.9968 

0.9987 

0.9966 

0.9966 

0.9985 

4.00 

0.9992 

0.9991 

0.9991 

0.9991 

0,9990 

0.9990 

4.10 

0.9994 

0.9994 

0.9994 

0.9994 

0.9993 

0.9993 

4.20 

0.9996 

0.9996 

0,9996 

0.9996 

0.9995 

0.9995 

4.30 

0.9997 

0.9997 

0.9997 

0.9997 

0,9997 

0.9997 

4.40 

0.9998 

0.9998 

0.9998 

0.9996 

0.9996 

0.9998 

4.50 

0.9999 

0.9999 

0.9999 

0.9999 

0.9999 

0,9999 

r 


-46- 


f 

13  TABLE  1 


c 

1.2000 

1.2500 

1.3000 

1.3500 

1.4000 

1.4500 

THETA 

Y 

2.0544 

2.1842 

2.3260 

2.4674 

2.6645 

2.8612 

0.0 

0.0400 

0.0359 

0.0311 

0.0265 

0.0222 

0.0162 

0.10 

0.0534 

0.0474 

0.0415 

0.0356 

0.0304 

0.0253 

0.20 

0.0668 

0.0616 

0.0545 

0.0476 

0.0409 

0.0346 

0.30 

0.0673 

0.0788 

0.0705 

0.0622 

0.0542 

0.0465 

0.40 

0.1002 

0.0995 

0.0898 

0.0601 

0.0706 

0.0614 

0.50 

0.1347 

0.1237 

0.1127 

0.1015 

0.0905 

0.0796 

0.60 

0.1630 

0.1518 

0.1394 

0.1268 

0.1142 

0.1017 

0.70 

0.1060 

0.1836 

0.1700 

0.1560 

0.1419 

0.1277 

0.80 

0.2335 

0.2192 

0.2045 

0.1893 

0.1737 

0.1579 

0.00 

0.2735 

0.2584 

0.2427 

0.2264 

0.2096 

0.1924 

1.00 

0.3164 

0.3007 

0.2643 

0.2671 

0.2493 

0.2309 

1.10 

0.3618 

0.3458 

0.3288 

0.3111 

0.2925 

0.2732 

1.20 

0.4001 

0.3929 

0.3756 

0.3577 

0.3387 

0.3167 

1.30 

0.4576 

0.4416 

0.4245 

0.4064 

0.3872 

0.3670 

1.40 

0.5066 

0.4909 

0.4742 

0.4564 

0.4374 

0.4172 

1.50 

0.5552 

0.5402 

0.5241 

0.5066 

0.4863 

0.4666 

1.60 

0.6030 

0.5888 

0.5735 

0.5570 

0.5393 

0.5203 

1.70 

0.6490 

0.6358 

0.6215 

0.6060 

0.5893 

0.5713 

1.60 

0.6028 

0.6807 

0.6675 

0.6532 

0.6377 

0.6210 

1.00 

0.7330 

0.7229 

0.7110 

0.6979 

0.6636 

0.6684 

2.00 

0.7710 

0.7621 

0.7514 

0.7397 

0.7270 

0.7131 

2.10 

0.8065 

0.7978 

0.7864 

0.7781 

0.7669 

0.7545 

2.20 

0.6375 

0.8301 

0.8219 

0.6129 

0.6031 

0.7923 

2.30 

0.6651 

0.8587 

0.6517 

0.6440 

0.6355 

0.8262 

2.40 

0.6601 

0.6837 

0.8778 

0.8713 

0.8642 

0.8563 

2.50 

0.0099 

0.9054 

0.9005 

0.6950 

0.8891 

Q .8825 

2.60 

0.9276 

0,9239 

0.9198 

0.9153 

0.9104 

0.9050 

2.70 

0.0424 

0.9394 

0.9361 

0.9325 

0.9285 

0.9240 

2.60 

0.9547 

0.9523 

0.9496 

0.9467 

0.9435 

0.9400 

2.00 

0.9647 

0.9628 

0.9608 

0.9584 

0.9559 

0.9531 

3.00 

0.9728 

0.9714 

0.9697 

0.9679 

0.9659 

0.9637 

3.10 

0.0703 

0.9762 

0.9769 

0.9755 

0.9740 

0.9723 

3.20 

0.0844 

0.9836 

0.9826 

0.9815 

0.9804 

0.9791 

3.30 

0.0664 

0.9877 

0.9870 

0.9862 

0.9853 

0.9844 

3.40 

0.9914 

0.9910 

0.9904 

0.9698 

0.9892 

0.9884 

3.50 

0.0936 

0.9934 

0.9930 

0.9926 

0.9921 

0.9916 

3.60 

0.9955 

0.9952 

0.9950 

0.9946 

0.9943 

0.9939 

3.70 

0.9068 

0.9966 

0.9964 

0.9962 

0.9959 

0.9956 

3.60 

0.9977 

0.9976 

0.9975 

0.9973 

0.9971 

0.9969 

3.00 

0.9984 

0.9983 

0.9982 

0.9961 

0.9980 

0.9978 

4.00 

0.9989 

0.9986 

0.9988 

0.9967 

0.9986 

0.9985 

4.10 

0.9993 

0.6992 

0.9992 

0.9991 

0.9990 

0,9990 

4.20 

0.9995 

0.9995 

0.9994 

0.9994 

0.9994 

0.9993 

4.30 

0.9997 

0.9996 

0.9996 

0.9996 

0.9996 

0.9995 

4.40 

0.9998 

0.9996 

0.9997 

0.9997 

0.9997 

0.9997 

4.50 

0.9999 

0.9998 

0.9998 

0.9998 

0.9996 

0.9998 
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TABLE  1 


c 

1.5000 

1.5500 

1.6000 

1.6500 

1.7000 

1.7500 

THETA 

V 

3.0602 

3.3243 

3.5966 

3.9011 

4.2419 

4.6240 

0.0 

0.0146 

0.0115 

0.0087 

0.0064 

0.0046 

0.0031 

0.10 

0.0207 

0.0165 

0.0129 

0.0097 

0.0071 

0.0050 

0.20 

0.0267 

0.0234 

0.0165 

0.0143 

0.0107 

0.0078 

0.30 

0.0302 

0.0324 

0.0262 

0.0207 

0.0159 

0.0118 

0.40 

0.0525 

0.0441 

0.0363 

0.0293 

0.0230 

0.0175 

0.50 

0.0601 

0.0589 

0.0494 

0.0405 

0.0325 

0.0254 

0.60 

0.0803 

0.0773 

0.0659 

0.0550 

0.0450 

0.0359 

0.70 

0.1136 

0.0997 

0.0662 

0.0732 

0.0610 

0.0498 

o.eo 

0.1421 

0.1263 

0.1107 

0.0956 

0.0611 

0.0675 

0.00 

0.1740 

0.1573 

0.1396 

0.1225 

0.1057 

0.0895 

1.00 

0.2120 

0.1928 

0.1735 

0.1541 

0.1350 

0.1164 

1.10 

0.2532 

0.2326 

0.2117 

0.1905 

0.1693 

0.1483 

1.20 

0.2070 

0.2764 

0.2541 

0.2314 

0.2064 

0.1853 

1.30 

0.3456 

0.323S 

0.3004 

0.2766 

0.2521 

0.2273 

1.40 

0.3059 

0.3735 

0.3499 

0.3254 

0.2999 

0.2738 

1.50 

0.4476 

0.4253 

0.4018 

0.3770 

0.3511 

0.3242 

1.60 

0.4909 

0.4782 

0.4551 

0.4306 

0.4046 

0.3777 

1.70 

0.5520 

0.5312 

0.5089 

0.4852 

0.4599 

0.4332 

1.80 

0.6026 

0.5833 

0.5623 

0.5397 

0,5154 

0.48^0 

1.00 

0.6516 

0.6337 

0.6142 

0.5930 

0.5702 

0.5457 

2.00 

0.6980 

0.6816 

0.6638 

0.6444 

0.6233 

0 . bO  06 

2.10 

0.7411 

0.7264 

0.7104 

0.6929 

0.6736 

0.6531 

2.20 

0.7805 

0.7676 

0.7534 

0.7379 

0.7209 

0.7024 

2.30 

0.8161 

0.8049 

0.7925 

0.7790 

0.7642 

0.7478 

2.40 

0.8476 

0.8380 

0.6275 

0.6159 

0.6031 

0.7890 

2.50 

0.8752 

0.6671 

0.8583 

0.6485 

0.637* 

0.6257 

2.60 

0.6989 

0.692J 

0.6849 

0.6766 

0.6677 

0.8577 

2.70 

0.9191 

0.9137 

0.9077 

0.9010 

0.6936 

0.8853 

2.60 

0.9360 

0.9316 

0.9266 

0.9214 

0.9i5u 

0.°087 

2.00 

0.9499 

0.9465 

0.9426 

0.9363 

0,9335 

0.9262 

3.00 

0.9613 

0.9566 

0.9555 

0.9522 

0.9484 

0.9442 

3.10 

0.9704 

0.9663 

0.9659 

0.9633 

0.9604 

0.9572 

3.20 

0.9776 

0.9760 

0,9742 

0.9722 

0.9700 

0.9675 

3.30 

0.9833 

0.9621 

0.9607 

0.9792 

0.9775 

0.975b 

3.40 

0.9876 

0.9867 

0.9657 

0.9646 

0.9833 

0.9819 

3.50 

0.9910 

0.9903 

0.9695 

0.9867 

0.9876 

0.9868 

3.60 

0.9935 

0.9930 

0.9924 

0.9916 

0.9912 

0.9904 

3.70 

0.9953 

0.9950 

0.9946 

0.9942 

0.9937 

0.9931 

3.90 

0.9V67 

0.9964 

0.9962 

0.9959 

0.9955 

0.9951 

3.00 

0.9977 

0.9975 

0.9973 

0.9971 

0.9969 

0.9966 

4.00 

0.9964 

0.9983 

0.9961 

0.9960 

0.9976 

0.9976 

4.10 

0.9989 

0,9968 

0.9987 

0.9986 

0.9965 

0.9984 

4.20 

0.9993 

0,9992 

0.9991 

0.9991 

0.9990 

0.9989 

4.30 

0.9995 

0.9995 

0.9994 

0.9994 

0.9993 

0.9993 

4.40 

0.9997 

0.9906 

0,9996 

0.9996 

0,9995 

0.9995 

4.50 

0.9998 

0.9998 

0,9997 

0.9997 

0.9997 

0.9997 
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35  TABLE  1 


c 

1 .8000 

1.8500 

1.9000 

1.9500 

2.0000 

2.0500 

THETA 

V 

5.0531 

5.5350 

6.0600 

6.6043 

7.3601 

8.1 764 

0.0 

0.0020 

0.0013 

0.0007 

0.0004 

0.0002 

0.0001 

0.10 

0.0034 

0.0021 

0.0013 

0.0007 

0.0004 

0.0002 

0.20 

0.0054 

0.0036 

0.0023 

0.0013 

0.0007 

0.0004 

0.30 

0 . 0084 

0.0058 

0.0038 

0.0024 

0.0014 

0.0007 

0.40 

0.0120 

0.0001 

0.0062 

0.0040 

0.0024 

0.0014 

O.SO 

0.0102 

0.0140 

0.0099 

0.0066 

0.0042 

0.0025 

0.60 

0.0270 

0.0210 

0.0152 

0.0106 

0.0070 

0.0044 

0.70 

0.0306 

0.0306 

0.0229 

0.0165 

0.0114 

0.0075 

o.eo 

0.0540 

0.0435 

0.0334 

0.0240 

0.0178 

0.0122 

0.00 

0.0743 

0.0603 

0.0476 

0.0365 

0.0270 

0.0192 

1.00 

0 .0085 

0.0817 

0.0661 

0.0521 

0.0308 

0.0293 

1.10 

0.1278 

0.1061 

0.0896 

0.0724 

0.0569 

0.0433 

1.20 

0.1624 

0.1400 

0.1184 

0.0961 

0.0792 

0.0622 

1.30 

0.2023 

0.1774 

0.1531 

0.1295 

0.1073 

0.0866 

1.40 

0.2471 

0.2202 

0.1934 

0.1670 

0.1415 

0.1173 

1.50 

0.2064 

0.2660 

0.2393 

0.2105 

0.1621 

0.1545 

1.60 

0.3404 

0.3201 

0.2000 

0.2594 

0.2287 

0.1982 

1.70 

0.4050 

0.3755 

0.3446 

0.3131 

0.2808 

0.2482 

1  •  BO 

0.4621 

0.4330 

0.4024 

0.3705 

0.3374 

0.3035 

1.00 

0.5105 

0.4915 

0.4617 

0.4302 

0.3972 

0.3629 

2.00 

0.5760 

0.5496 

0.5212 

0.4910 

0.4589 

0.4251 

2.10 

0.6305 

0.6061 

0.5797 

0.5513 

0.5209 

0.4884 

2.20 

0.6821 

0.6600 

0.6359 

0.6006 

0.5816 

0.5512 

2.30 

0.7200 

0.7103 

0.6666 

0.6653 

0.6307 

0.6119 

2.40 

0.7735 

0.7564 

0,7375 

0.7166 

0.6941 

0.6693 

2.50 

0.6124 

0.7978 

0.7616 

0.7637 

0.7440 

0.7222 

2.60 

0.8466 

0.6343 

0.6206 

0.6055 

0,7686 

0.7700 

2.70 

0.6762 

0.6660 

0.6546 

0.6420 

0.6279 

0.6123 

2.60 

0.0013 

0.6930 

0.6637 

0.6734 

0.6618 

0.6469 

2.00 

0.0222 

0.0156 

0.0062 

0.6996 

0.6905 

0.6801 

3.00 

0.0305 

0.0343 

0.0264 

0.9216 

0.9144 

0.9061 

3.10 

0.0535 

0.0494 

0.0446 

0.9397 

0.9330 

0.9274 

3.20 

0.0647 

0.0616 

0,0560 

0.9541 

0.0406 

0.9446 

3.30 

0.0735 

0.0711 

0.0665 

0.0655 

0.0621 

0.9582 

3.40 

0.0803 

0.0786 

0.0766 

0.0743 

0.0718 

0.0689 

3.50 

0.0856 

0.0643 

0.0628 

0.9811 

0.9793 

0.9772 

3.60 

0.9606 

0.9866 

0.0675 

0.9663 

0.9849 

0.9834 

3.70 

0.9025 

0.0918 

0.091 1 

0.9902 

0.9892 

0.0861 

3. BO 

0.0047 

0.0942 

0.0937 

0.0930 

0.9923 

0.9915 

3.00 

0.9963 

0.0050 

0.0955 

0.0951 

0.9946 

0.9941 

4.00 

0.0974 

0.9072 

0.0960 

0.9966 

0.9963 

0.0959 

4.10 

0.0062 

0.0061 

0.9979 

0.9977 

0.9974 

0.9972 

4.20 

0.0968 

0.9967 

0.9966 

0.0964 

0.9983 

0.9981 

4.30 

0.0002 

0.9991 

0.0990 

0.9960 

0.9986 

0.0987 

4.40 

0.0005 

0.0094 

0,9994 

0.9993 

0.0992 

0.9991 

4.50 

0.0007 

0.0096 

0.0996 

0.0995 

0.0995 

0.9994 
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TABLE  1 


c 

2.1000 

2.1500 

2.2000 

2.2500 

2.3000 

2.3500 

THETA 

V 

0.0703 

10.0870 

11.2459 

12.5692 

14.0634 

15.8196 

0.0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.10 

0.0001 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.20 

0.0002 

0.0001 

0.0000 

0.0000 

0.0000 

0.0000 

0.30 

0.0004 

0.0002 

0.000! 

0.0000 

0.0000 

0.0000 

0.40 

0.0007 

0.0004 

0.0002 

0.0001 

0.0000 

0.0000 

o.so 

0.0014 

0.0007 

0.0003 

0.0001 

0.0001 

0.0000 

0.60 

0.0026 

0.0014 

0.0007 

0.0003 

0.0001 

0.0000 

0.70 

0.0046 

0.0027 

0.0014 

0.0007 

0.0003 

0.0001 

o.eo 

0 . 00  7e 

0.0048 

0.0028 

0.0015 

0.0007 

0.0003 

0.00 

0.0130 

0.0064 

0.0051 

0.0029 

0.0015 

0.0007 

1.00 

0.0207 

0.0140 

0.0089 

0.0053 

0.0030 

0.0015 

1.10 

0.0318 

0.0224 

0.0150 

0.0095 

0.0056 

0.0031 

1.20 

0.0473 

0.0346 

0.0242 

0.0161 

0.0101 

0.0059 

1.30 

0.0680 

0.0516 

0.0376 

0.0263 

0.0174 

0.0109 

l.ao 

0.0947 

0.0743 

0.0564 

0.0411 

0.0286 

0.0189 

1.50 

0.1282 

0.1037 

0.0614 

0.0618 

0.0450 

0.0313 

1.60 

0.1686 

0.1402 

0.1136 

0.0694 

0.0679 

0.0495 

1.70 

0.2157 

0.1839  • 

0.1534 

0.1246 

0.0983 

0.0749 

1.80 

0.2690 

0.2346 

0.2007 

0.1680 

0.1370 

0.1064 

1.00 

0.3275 

0.2914 

0.2551 

0.2191 

0.1842 

0.1509 

2.00 

0.3897 

0.3530 

0.3154 

0.2773 

0.2393 

0.2021 

2.10 

0.4540 

0.4177 

0.3800 

0.3410 

0.3013 

0.2614 

2.20 

0.5186 

0.4838 

0.4470 

0.4064 

0.3683 

0.3272 

2.30 

0.5816 

0.5493 

0.5145 

0.4775 

0.4384 

0.3974 

2.40 

0.6421 

0.6126 

0.5806 

0.5460 

0.5090 

0.4697 

2.50 

0.6983 

0.6720 

0.6433 

0.6120 

0.5781 

0.5415 

2.60 

0.7494 

0.7266 

0.7015 

0.6739 

0.6436 

0.6106 

2.70 

0.7949 

0.7755 

0.7541 

0.7303 

0.7040 

0.6749 

2.80 

0.8345 

0.6184 

0.8005 

0.7604 

0.7581 

0.7333 

2.00 

0.8664 

0.8552 

0.6405 

0,8240 

0.8055 

0.7846 

3.00 

0.8967 

0.6862 

0.6744 

0.6610 

0.8460 

0.6292 

3.10 

0.9200 

0.9118 

0.9024 

0.8919 

0.8799 

0.8664 

3.20 

0.9389 

0.9325 

0.9252 

0.9170 

0.9077 

0.6971 

3.30 

0.9539 

0.9490 

0.9435 

0.9372 

0.9300 

0.9219 

3.40 

0.9657 

0,9620 

0.9578 

0.9531 

0,9477 

0.9415 

3.50 

0.9748 

0.9720 

0.9689 

0.9654 

0.9614 

0.9568 

3.60 

0.9816 

0.9797 

0.9774 

0.9746 

0.9719 

0.9685 

3.70 

0.9668 

0.9854 

0.9837 

0,9619 

0.9797 

0.9773 

3.80 

0.9906 

0.9896 

0.9884 

0.9671 

0.9856 

0.9839 

3.00 

0.9934 

0.9927 

0.9919 

0.9909 

0.9699 

0.9887 

«.oo 

0.9954 

0.9949 

0.9944 

0.9937 

0.9930 

0.9921 

4.10 

0.9969 

0.9965 

0.9961 

0,9957 

0.9952 

0.9946 

4.20 

0.9979 

0.9976 

0,9974 

0.9971 

0.9967 

0.9963 

4.30 

0.9986 

0.9984 

0.9982 

0.9960 

0.9978 

0.9975 

4.40 

0.9991 

0.9989 

0.9966 

0.9987 

0.9985 

0.9984 

4.50 

0.9994 

0.9993 

0.9992 

0.9991 

0.9990 

0.9989 
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1 .  INTRODUCTION 

It  has  been  observed  empirically  (Gringorten, 
1966,  1968)  that  fluctuations  in  upper-level  wind 
velocity,  rainfall,  temperature,  and  other  meteor¬ 
ological  variates  may  be  modeled  by  a  type  of  ran¬ 
dom  Markov  process.  A  good  fit  is  obtained  over 
a  broad  range  of  values  when  the  so-called  Omstein- 
Uhlenbeck  process  described  below  is  employed. 
Acceptance  of  the  validity  of  such  a  model  permits 
utilization  of  its  theoretical  predictions  for 
long-range  operational  planning,  with  on-line  pre¬ 
dictions  available  from  software.  Such  software 
is  in  development  at  the  U.S.A.F.  Geophysical 
Research  Laboratory. 

The  O-U  process  is  the  simplest  random  pro¬ 
cess  having  a  stationary  distribution  whose  sample 
paths  are  continuous  functions  of  time.  It  is 
correspondingly  easiest  to  work  with.  Certain 
passage  time  distributions  needed  for  meteoro¬ 
logical  planning  have  been  difficult,  however,  to 
calculate,  and  have  only  recently  become  available 
(Keilson  and  Ross,  1975). 

The  random  variable  of  central  interest  to 
this  paper  is  M(8) ,  the  maximum  of  an  0-U  vaTiate 
over  an  interval  of  duration  6  in  the  absence  of 
any  conditioning  history.  The  distribution  of 
the  maximum  M(e)  relates  directly  to  the  passage 
time  distributions  in  the  Keilson-Ross  paper. 

The  calculation  technique  employed  there  evaluated 
zeros  and  residues  in  the  complex  plane  of  the 
parabolic  cylinder  functions,  and  such  techniques 
again  feature  in  tabulating  the  distribution  of 
M(8). 

This  paper  describes  briefly  the  results 
obtained  for  the  interval  maximum.  A  more  com¬ 
plete  and  detailed  presentation  is  available 
(Keilson  and  Ross,  1978). 

It  must  be  emphasized  that  the  0-U  model  as¬ 
sumes  the  negligibility  of  diurnal  and  seasonal 
variation.  Where  such  influence  cannot  be  neg¬ 
lected,  the  results  described  here  cannot  be 
applied  without  modification. 

2.  THE  0-U  PROCESS 

The  Omstein-Uhlenbeck  process  (Cox  and 
Miller,  1965  and  Keilson  and  Ross,  1975)  describes 
wide  band  fluctuations  averaged  by  exponential 
smoothing.  For  the  resulting  Markov  process,  X(t), 
starting  at  X(0)  *  Xq,  one  has  in  the  absence  of 

observation  subsequent  to  t  «  0, 


H.  F.  Ross 

The  Medical  School,  University  of  Birmingham 
Birmingham,  England  B15  2TJ 


E[X(t))  -  x0e_t  ,  t  >  0  ;  (1) 

Var[X(t)]  -  1  -  e'2t  ,  t  >  0  .  (2) 

As  time  progresses,  X(t)  settles  into  a  stationary 
distribution  which  is  standard  normal,  i.e.,  one 
has 

2 

x  '  f- 

lim  P[X(t)  s  x]  -  /  ^ - dy  .  (J) 

t->“  -•  /TH 

Of  central  interest  is  the  stationary  0-U 
process  Xs(t)  in  which  the  variate  X(t)  starts  in 

its  stationary  normal  distribution  and  keeps  tf'is 
distribution  for  all  time.  The  maximum  M(B)  is 
defined  in  terms  of  this  stationary  process  by 

M(e)  *  max  X  (t ' )  ,  any  t  >  0  .  (4) 

tSt’St*B  5 

Its  distribution  will  be  designated  by 
F0(y)  ■  P[M(6)  s  y]  .  P[Xs(f)  S  y; 

t  s  t’  s  t  ♦  e)  . 

The  stationary  0-U  process  is  Gaussian,  and 
is  sometimes  called  the  Gaussian  Markov  process. 

It  has  for  its  covariance,  when  tj  <  t^, 

-(t  -t  ) 

cov[X(tj),  X(t2)]  -  e  1  .  (5) 

All  correlation  is  therefore  lost  when 
tj  -  tj  »  1.  The  process  has  a  relaxation  time 

or  forgetting  time  (Keilson,  1979) 

T  •  1  .  (6) 

3.  THE  MODELING  OF  METEOROLOGICAL  VARIATES 

A  meteorological  process  Z(t)  siodeled  via  an 
0-U  process  X(t)  has  three  parameters: 

(a)  the  mean  value  m  of  the  associated  stationary 
process  Zg (t ) ; 

(b)  the  standard  deviation  o  of  the  associated 
stationary  process  Zs(t); 
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(c)  the  relaxation  tine  T  of  the  associated  sta¬ 
tionary  process  Zs(t). 

The  process  Z(t)  then  relates  to  the  under¬ 
lying  0-U  process  X(t)  by 

Z(t)  «  m  ♦  o  X(£)  .  (7) 

These  three  paraneters  will  vary  with  geo¬ 
graphical  location  and  altitude,  apart  from  diur¬ 
nal  and  seasonal  changes.  The  use  of  these  three 
parameters  is  required  for  application  of  the 
graphs  and  tables.*  How  they  are  used  will  be 
discussed  when  the  graphs  are  presented. 

4.  INADEQUACY  OF  EXTREME  VALUE  DISTRIBUTIONS 

The  extreme  value  distributions  of  Gnedenko, 
Gumbel  and  others  are  known  to  be  valid  for  M(6), 
our  interval  maximum  for  the  0-U  process,  when  6 
is  sufficiently  large,  as  shown,  for  example,  by 
Darling  and  Erdos,  1955.  The  validity  of  such 
asymptotic  description  has  been  known  to  set  in 
very  slowly,  and  the  error  in  the  asymptotic  dis¬ 
tribution  of  M(8)  for  any  value  of  8  has  not  been 
known.  The  results  presented  here  are  accurate 
for  the  0-U  model.  They  permit  comparison  with 
extreme  value  prediction  and  confirm  the  inade¬ 
quacy  of  asymptotic  extreme  value  theory  for  most 
interval  durations  8  of  operational  planning 
interest . 

5.  NUMERICAL  RESULTS  AND  GRAPHS 

The  distribution  of  the  maximum  M(6)  is 
directly  related  to  the  passage  time  distribution 
for  the  process  Xs(t).  This  arises  from  the  ob¬ 
servation  that  the  event  that  the  maximum  value 
of  a  process  sample  is  less  than  y  in  an  interval 
of  duration  8  is  the  same  as  the  event  that  the 
process  sample  takes  an  interval  longer  than  8  to 
exceed  a  value  y.  Consequently,  if  the  passage 
time  distribution  of  the  process  is  denoted  by 


This  form  of  presentation  of  the  graph  is  most 
useful  if  the  type  of  question  being  asked  is 
"given  one  of  the  threshold  values  y  «  0.5,  1.0, 
etc.,  what  is  the  probability  of  exceeding  that 
threshold  in  a  time  interval  of  given  length?". 

Though  the  graphs  terminate  at  50  relaxation 
times,  they  may  easily  be  extended  further  since 
beyond  that  point  they  are  adequately  approxi¬ 
mated  by  exponential  curves,  which  means  that 
values  of  Sy(6)  for  8  >  50  may  be  obtained  from 

the  equation  Sy(6)  «  1  -  e  ,l^)e 

Table  I  gives  the  value  of  X(y)  for  the  set 
of  values  of  the  y  parameter  used  in  the  graph. 


Table  I 


y 

X(y) 

0.5 

0.6488 

1.0 

0.3882 

1.5 

0.2087 

2.0 

0.09727 

2.5 

0.03768 

3.0 

0.01161 

3.5 

0.002752 

Fig.  2  shows  F0(y)  as  a  function  of  y  for 

various  values  of  c,  where  c  ■  ^21og6.  This  form 
of  presentation  is  most  appropriate  for  showing 
the  way  in  which  the  y  and  8  parameters  interact 
for  large  value  of  8  (long  time  durations) .  The 
independent  variable  c  *  /21og6  has  been  used 
because  it  can  be  shown  that  in  the  limit,  for 
large  values  of  8,  the  variable  c  is  a  good  meas¬ 
ure  of  location  for  the  distribution  F.(y).  This 

feature  can  be  seen  in  the  graph  even  for  the  rel¬ 
atively  small  values  of  c  used.  The  relation 
between  c  and  8  is  given  explicitly  in  Table  II 
for  the  values  of  c  used  in  the  graph. 


S  (8)  *  P[the  time  taken  for  the  process  to  first 
reach  a  value  y  is  s  8] 

then 

1  -  S  (8)  *  P[that  the  process  does  not  exceed  a 
'  value  y  within  a  time  interval  of 

duration  8] 

'  Fe(>}  • 

Fig.  1  shows  Sy(8)  as  a  function  of  8  for  various 
values  of  y.  The  time  units  are  relaxation  times 
of  the  process  under  consideration  in  accordance 
with  the  discussion  in  Section  S.  Note  that  at 
the  left  edge  of  the  graph  the  curves  do  not  go 
to  0  but  to  values  given  by  the  expression 

2 

Sy(0)  -  -L  5  e'  2"  -  •  («) 

y  Ji 7  — 


•Tables  are  given  in  Keilson  and  Ross,  1978. 


Table  II 


c 

o 

t> 

N 

a> 

1.0 

1.65 

1.5 

3.08 

2.0 

7.39 

2.5 

2.28  x 

3.0 

9.00  x 

3.5 

4.57  x 

4.0 

2.98  x 

4.5 

2.49  x 

5.0 

2.68  x 

0.5c 


10 

10 

102 

103 

104 

106 


For  large  values  of  e.  F  Cy)  approaches  the 
8 

extreme  value  distribution  as  a  limiting  form. 

This  approach  is  rather  slow,  however.  Numerical 
calculations  indicate  that  it  is  not  until  c  reaches 
a  value  of  about  5  that  the  error  in  the  distribu¬ 
tion  function  from  using  the  extreme  value  form  is 
down  to  approximately  1%.  Table  II  shows  that  that 
value  of  c  corresponds  to  2.68  x  103  relaxation 
times,  which  for  a  relaxation  time  of  one  now  means 
a  period  of  about  30  years. 
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